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1.  This  is  the  Final  Report  on  the  contract.  It  describes  the  analyti- 
cal work  performed  to  modify  the  RADC  Polar  Ionospheric  Model  on  the 
basis  of  measurements  made  with  the  DNA  Project  617  Incoherent -scatter 
radar  at  Chatanika,  Alaska.  This  radar  had  the  unique  capability  of 
measuring  the  electron  density  profile  in  the  upper  E-reglon  when  a 
local  minimum  or  "valley"  is  present. 

2.  The  objective  of  this  work  was  an  Improved  specification  of  the 
auroral  E-layer  and  the  altitude  interval  between  the  E and  F regions. 

3.  This  work  is  valuable  because  of  its  application  to  ducted  propaga- 
tion studies.  The  altitude  interval  between  the  E and  F regions  some- 
times acts  as  a duct  in  which  High  Frequency  signals  may  travel  for 
great  distances  without  returning  to  earth.  Ducted  rays  are  less 
attenuated  than  normal  modes  because  they  make  fewer  traversals  of  the 
D region,  where  most  of  the  absorption  occurs. 

4.  The  results  will  be  used  in  ray-tracing  calculations  to  determine 
probability  of  ducting  for  selected  paths  in  the  design  of  long  range 
propagation  experiments. 


WILLIAM  F.  RING  / 

Contract  Monitor 
Propagation  Branch 
Electromagnetic  Sciences  Division 
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SUMMARY 


The  ability  to  predict  the  propagation  of  HF  radio  waves  along 
transauroral  paths  requires  a model  of  the  spatial  distribution  of 
electron  density  in  the  hlgh-latltude  ionosphere.  An  empirical  mean 
model  has  been  developed  by  Elkins  and  coworkers  (Elkins  and  Rush,  1973) 
and  is  referred  to  as  the  RADC-POLAR  model.  This  model  has  been  used  by 
RADC  and  others  for  predictive  studies  with  three-dimensional  raytracing 
programs . 

An  ionospheric  region  of  particular  Interest  for  transauroral 
propagation  is  the  altitude  Interval  between  the  E and  F regions.  This 
region  sometimes  acts  as  a duct  in  which  HF  signals  may  travel  for  great 
distances  without  returning  to  earth.  The  presence  of  a "valley",  or 
local  minimum  above  the  E layer,  is  favorable  (but  not  essential)  for 
ducted  propagation.  Rays  that  are  ducted  will  be  less  attenuated  than 
normal  modes  because  they  make  fewer  traversals  of  the  D region,  where 
most  of  the  absorption  occurs.  Such  propagation  modes  are  especially 
important  during  magnetically  disturbed  periods,  when  auroral  absorption 
can  completely  obliterate  an  HF  signal. 

The  RADC-POLAR  model  is  based  primarily  on  ground-based  and 
satellite-borne  vertical-incidence  ionosonde  data  that  do  not  measure 
the  density  profile  in  the  upper  E-reglon  when  a valley  is  present. 
However,  complete  profiles  of  electron  density,  from  the  lower  E region 
to  well  above  the  F-reglon  peak,  can  be  measured  by  an  incoherent -scat ter 
radar.  The  purpose  of  our  study  was  to  analyze  measurements  by  the  DNA 
Project  617  Incoherent -scatter  radar  at  Chatanika,  Alaska,  and  to  extract 
a synoptic  model  that  describes  the  characteristic  variations  of 
ionization  at  auroral  latitudes. 
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Because  the  existing  RADC  model  was  derived  from  a larger  and 
more  geographically  extensive  data  base  than  that  used  in  this  study, 
many  of  its  features  have  been  retained  in  the  new  model.  The  major 
modification  that  we  sought  was  an  improved  specification  of  the  auroral 
E layer  and  the  valley  above  it. 

The  data  used  in  this  work  consisted  of  21  sets,  each  of  24-hours 
duration,  and  comprising  approximately  5,000  profiles.  Data  were 
selected  to  give  a representative  distribution  with  season  and  geo- 
magnetic activity.  Two  types  of  data  were  used:  elevation  scans  in  the 
magnetic  meridian  plane,  and  near-vertical  measurements  over  Chatanika . 
When  the  radar  is  operated  in  an  elevation-scan  mode,  it  measures  the 
distribution  of  ionization  in  both  altitude  and  latitude  with  a spatial 
resolution  of  10  km  over  a horizontal  extent  of  more  than  600  km.  Two 
days  of  such  data  were  examined  in  order  to  obtain  a qualitative  descrip- 
tion of  the  diurnal  morphology  of  the  high-latitude  ionosphere.  The 
remainder  of  the  data  consisted  of  electron  density  profiles  in  the 
region  nearly  overhead  of  Chatanika. 

An  examination  of  the  elevation-scan  data  showed  that,  as  expected, 
the  ionization  could  be  attributed  to  two  sources:  solar  photoionization 
md  particle  precipitation.  The  dayside  F-reglon  profile  appears  to  be 
consistent  with  production  by  solar  ionization.  However,  the  daytime 
E-reglon  profile  is  often  observed  to  have  a dense  E layer  that  is 
evidently  due  to  particle  precipitation  south  of  the  dayside  auroral 
oval.  A',  night  the  auroral  E-layer  has  spatial  and  temporal  variations 
characteristic  of  diffuse  and  discrete  auroral  precipitation.  Often  the 
electron  density  decreases  monotonlcally  with  altitude  above  the  E layer 
with  no  indication  of  an  F layer  or  valley. 

From  each  of  the  electron  density  profiles  we  extracted  the 
fundamental  parameters  used  in  the  RADC-POMR  model  to  derive  the  model 
profiles.  Among  these  parameters  are  the  height  and  maximum  density  of 
the  auroral-E  layer,  the  height  and  minimum  density  in  the  valley  (if 
present)  above  the  E-reglon  peak,  and  the  height,  density,  and  semi- 
thickness  of  the  FI  and  F2  layers.  A statistical  analysis  of  these 
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parameters  wts  made  to  determine  their  characteristic  variation  with 
solar  zenith  angle  and  magnetic  activity.  Solar  zenith  angle  was 
chosen  as  a sorting  parameter  In  order  to  combine  diurnal  and  seasonal 
variations  Into  a single  index. 

The  results  of  this  statistical  analysis  were  used  to  modify  the 

RADC-POIAR  model.  The  modifications  were  made  so  as  to  be  compatible 

with  the  present  predictive  model,  which  defines  auroral  features  as 

functions  of  the  planetary  geomagnetic  index  K . We  found  that  the 

P 

profile  parameters  were  only  weakly  correlated  with  and  that  a 

significantly  stronger  correlation  exists  with  the  local  magnetic  index, 

K.  Since  we  wanted  to  retain  K as  the  descriptor  of  magnetic  activity, 

P 

a relationship  was  sought  between  the  planetary  index  and  the  local  K 
values  at  various  latitudes  and  local  times.  A meaningful  relationship 
was  derived  from  geomagnetic  data  obtained  over  an  entire  year  at  16 
high-latitude  stations.  Incorporation  of  this  relationship  into  the  model 
Increased  the  spatial  and  temporal  structure  within  the  auroral  oval. 

Other  modifications  that  were  made  to  the  RADC-POLAR  model  include  the 
addition  of  four  new  profile  shapes,  an  increase  in  the  size  of  the 

nightslde  oval,  the  provision  for  particle-produced  ionization  on  the  i 

dayside  south  of  the  auroral  oval,  and  a reduction  in  the  electron  i 

density  of  the  midlatitude  auroral  E-layer.  j 

Our  analysis  is  similar  to  that  used  in  the  development  of  the  j 

RADC-POLAR  model,  in  that  both  produce  statistical  mean  models.  We  did  j 

not  model  small-scale  (<100  km)  structures  or  phenomena  that  were  | 

transient  on  a time-scale  smaller  than  the  3-hour  K index.  As  a result  ! 

P j 

the  features  in  our  model  are  characteristic  of  persistent,  widespread  j 

phenomena,  such  as  the  diffuse  aurora  and  the  dayside  "drizzle"  zones. 

Features  such  as  auroral  arcs,  the  diffuse  aurora  "boundary  feature," 

and  auroral  substorms  are  not  included.  However,  the  use  of  a local  K 

that  varies  with  latitude  and  local  time  for  a specified  did  result 

in  an  auroral  ionospheric  model  with  more  structure  than  that  in  the 

original  model  and,  presumably,  a more  realistic  one. 


i 
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To  demonstrate  the  compatibility  of  the  new  model  with  the  RADC 
raytracing  code,  a number  of  rays  were  traced  from  a point  in  the  north- 
western U.S.  toward  Europe.  Rays  were  launched  at  the  same  angles 
through  the  unmodified  ionospheric  model  to  Illustrate  the  effect  of 
differences  between  the  two  versions. 

In  addition  to  the  model  of  the  refracting  regions,  a model  of 
auroral  absorption  was  developed  for  possible  incorporation  into  the 
RADC  code  at  a later  time.  This  model,  developed  from  previously  compiled 
riometer  datg,  describes  the  variation  of  median  auroral  absorption  as  a 
function  of  corrected  geomagnetic  latitude,  longitude,  season,  local 
time,  solar  activity,  and  magnetic  activity. 

The  primary  end  product  of  our  work  is  the  modified  ionospheric 
model.  The  major  components  of  it  are  listed  in  Appendix  C,  and  a 
complete  operating  program  is  being  sent  to  RADC  separate  from  this 
report . 


1.  INTRODUCTION 


The  ability  to  predict  the  propagation  of  HF  radio  waves  along 
transauroral  paths  requires  a model  of  the  spatial  distribution  of 
electron  density  In  tije  hlgh-latltude  Ionosphere.  For  this  purpose, 
Rlklns  { nd  coworkers  'Elkins  and  Rush,  1973)  developed  a model  referred 
to  as  th j RADC-1976-P(JLAR  Ionosphere.  This  report  describes  modifica- 
tions to  that  model  recommended  by  SRI  on  the  basis  of  measurements  made 
with  the  Incoherent-scatter  radar  at  Chatanlka,  Alaska. 

An  Ionospheric  feature  of  particular  Interest  for  transauroral 
propagation  Is  the  altitude  Interval  between  the  E and  F regions.  This 
region  sometimes  acts  as  a duct  In  which  HF  signals  may  travel  for  great 
distances  without  returning  to  earth.  Ducted  signals  will  be  relatively 
unattenuated  because  they  make  fewer  traversals  of  the  D i-eglon,  where 
most  of  the  absorption  occurs.  Such  modes  are  especially  important 
during  magnetically  disturbed  periods,  when  auroral  absorption  can 
completely  obliterate  an  HF  signal  as  it  traverses  the  disturbed  D 
region. 

The  RADC-POLAR  model  is  based  primarily  on  lonosonde  data  (both 
ground-based  and  satellite-borne)  that  do  not  measure  the  density 
profile  in  the  upper  E region  when  a local  minimum,  or  "valley",  is 
present.  In  order  to  determine  accurately  the  electron  density  there. 

It  Is  necessary  to  use  rocket-borne  Instrumentation  or  an  incoherent- 
scatter  radar.  The  purpose  of  this  study  is  to  analyze  measurements  by 
the  incoherent-scatter  radar  at  Chatanika , Alaska,  so  as  to  describe  the 
characteristic  variation  of  electron  density  at  auroral  latitudes. 

The  data  used  In  this  study  consisted  of  21  sets,  each  of  24-hours 
duration,  and  comprising  approximately  5,000  electron  density  profiles. 
Data  were  selected  to  give  a representative  distribution  of  season  and 
magnetic  activity.  Two  types  of  data  were  used.  Data  in  an  elevation 


scan  mode  provided  Information  about  the  distribution  of  electron  density 
with  both  altitude  and  latitude.  These  data  were  examined  in  order  to 
obtain  a qualitative  description  of  the  morphology  of  the  high-latitude 
ionosphere.  The  remainder  of  the  data  (approximately  90%  of  the  total) 
were  measurements  of  electron  density  profiles  restricted  to  the  region 
nearly  overhead  of  Chatanika. 

From  each  of  the  electron  density  profiles  we  extracted  the 
parameters  required  in  the  RADC-POLAR  model  to  specify  the  high- latitude 
Ionosphere.  A statistical  analysis  of  these  parameters  was  made  to 
determine  their  characteristic  variation  with  solar  zenith  angle  and 
magnetic  activity.  The  results  of  this  analysis  were  used  to  modify 
the  RADC-POLAR  model. 

As  demonstrations  of  the  use  of  the  modified  ionospheric  model, 
raytraclug  calculations  were  made  along  a particular  propagation  path. 
Raytracings  were  also  made  through  the  unmodified  ionospheric  model  as 
lllustra blons  of  the  effect  of  differences  between  the  two  versions. 

In  addition  to  the  model  of  ionospheric  electron  density,  we  also 
developed  a model  of  auroral  absorption.  Previously  complied  riometer 
data  were  used  to  synthesize  a description  of  the  variation  of  auroral 
absorption  with  latitude,  local  time,  and  magnetic  activity. 

The  report  is  organized  in  the  following  way.  Section  2 describes 
the  Chatanika  radar  and  its  measurements  of  electron  density.  Then  we 
specify  the  data  base  used  in  this  study  and  the  rationale  for  its 
selection.  In  Section  3 representative  examples  of  these  data  are  shown, 
together  with  a description  of  the  morphology  of  the  high-latitude 
ionosphere.  The  extraction  of  profile  parameters  and  their  statistical 
variations  are  presented  in  Section  4.  There  we  also  explain  the 
modifications  that  were  made  to  the  RADC-POLAR  model,  and  then  compare 
the  predictions  of  the  models  to  actual  Chatanika  data.  Examples  of 
raytracings  through  the  model  ionospheres  are  presented  in  Section  5. 

The  auroral  absorption  model  is  described  in  Section  6 and  the  conclusions 


of  our  study  are  given  in  Section  7,  Appendices  A through  C describe 
temperature  corrections  for  Chatanika  electron  density  measurements,  the 
relationship  between  and  local  K,  and  the  actual  modifications  to  the 
RADC-POLAR  software. 


2.  IONOSPHERIC  MEASUREMENTS  WITH  THE  CHATANIKA  RADAR 


The  electron  density  profiles  used  in  this  study  are  based  on 
measurements  made  with  the  Defense  Nuclear  Agency  Project  617  incoherent- 
scatter  radar  at  Chatanika,  Alaska.  This  radar,  generally  referred  to  as 
simply  the  Chatanika  radar,  has  been  operated  in  Alaska  by  SRI  since 
June  1971  and  is  the  only  incoherent-scatter  radar  at  a high  latitude. 

The  geographical  location  of  the  radar  is  65.1°N,  145. 5°W,  corresponding 
to  a geonagnetlc  location  of  64. 9° N,  103. 3° W and  an  L of  5.7. 

The  radar  transmits  a high-power  (typically  3.5  MW)  pulse  of  UHF 
(1290  MHz)  electromagnetic  radiation  and  the  backscattered  signal  from 
the  ionosphere  is  analyzed  for  both  amplitude  and  frequency  distribution. 
This  information  is  then  used  to  compute  the  electron  density,  electron 
and  ion  temperature,  plasma  drift,  and  other  derivable  parameters  such 
as  the  ionospheric  electric  field,  conductivity,  current,  and  neutral 
wind.  The  Chatanika  radar  and  associated  data -analysis  techniques  have 
been  described  in  a number  of  publications  (Leadabrand  et  al.,  1972; 

Baron,  1974,  1977;  Wickwar , 1975).  In  this  report  we  will  describe  only 
the  techniques  used  to  measure  the  spatial  variation  of  electron  density 
and  the  Chatanika  data  set  used  In  this  study. 

2.1  Electron  Density  Measurements 

Incoherent  scatter  is  basically  a technique  for  measuring  ionospheric 
electron  concentrations  and  other  parameters  by  using  a ground-based 
radar  facility.  Electrons  in  the  ionosphere,  when  excited  by  a UHF  radio 
wave,  will  backscatter  a small  amount  of  energy.  The  amount  of  energy 
scattered  is  proportional  to  the  radar  cross  section  of  a single  electron 
and  the  number  of  electrons  present  in  the  scattering  volume  (defined  by 
the  rada  • beam  and  the  transmitted  pulse  length) . By  measuring  the  back- 
scattered  power,  one  is  able  to  calculate  the  electron  number  density. 
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The  electron  density  is  normally  computed  by  assuming  that  the 
electron  cross  section  is  equal  to  the  classical  Thomson  cross  section. 

I This  result  is  generally  not  exact  because  the  radar  cross  section  of 

I the  electrons  is  not  a constant;  due  to  plasma  interactions  it  is  a 

function  of  the  electron  and  ion  kinetic  temperatures,  and  to  a smaller 
extent,  a function  of  the  ratio  of  the  Debye  length  to  the  operating 
wavelength  (Evans,  1969),  The  effect  of  neglecting  this  correction  factor 
that  accounts  for  plasma  interactions  is  to  always  underestimate  the  true 
electron  density.  This  correction  factor  is  generally  negligible  in  the 

j E-reglon  (<  ]50  km  altitude).  At  higher  altitudes  the  correction  often 

I 

I can  be  a factor  of  2 or  more.  However,  because  the  correction  factor  is 

generally  slowly  varying  with  altitude  there  is  only  a small  distortion 
in  the  shape  of  the  electron  density  profile  in  the  upper  E and  lower  F 
region.  The  effect  of  the  plasma  correction  factor  is  discussed  in  more 

j detail  in  Appendix  A. 

Measurements  of  E-  and  F-region  electron  densities  at  Chatanlka 
are  normally  done  using  short  (40  to  67  |is)  transmitted  pulses  with  192 
receiver  range  gates.  The  pulse  length  results  in  a spatial  range 
resolution  of  6 to  10  km.  These  data  are  processed  with  signal 
Integration  periods  of  1 s or  more,  the  choice  being  dependent  on  the 
magnitude  of  the  electron  densities  and  the  dynamic  character  of  the 
medium. 

Because  of  its  altitude  resolution  of  6 to  10  km,  the  Chatanlka 
radar  is  unable  to  detect  structures  of  smaller  scale  size  that  may  be 
responsible  for  sporadic  E,  Although  sporadic  E is  important  to  HF 
propagation,  we  made  no  attempt  to  include  it  in  the  model. 

The  Chatanlka  radar  has  a high-gain  antenna  with  a 3-dB  beamwidth 
of  0.5°,  As  a result,  when  the  antenna  is  held  stationary  the  horizontal 
spatial  resolution  in  the  E region  is  approximately  1 km.  Because  the 
antenna  Is  fully  steerable,  a variety  of  antenna  opei'ating  modes  are 
employed.  In  order  to  obtain  vector  estimates  of  plasma  drift,  measure- 
ments are  made  at  several  pointing  directions.  Generally,  the  radar 
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is  pointed  at  some  elevation,  typically  60°  to  70° , and  then  is  either 
scanned  in  azimuth  or  sequentially  pointed  to  three  separate  azimuths. 
About  six  minutes  are  required  for  each  cycle.  For  the  purpose  of  this 
study,  measurements  in  this  mode  are  referred  to  as  overhead  data 
because  the  horizontal  separation  from  Chatanlka  is  less  than  100  km. 

An  alternative  mode  of  operation  is  to  scan  the  radar  mechanically 
in  elevation  with  the  azimuth  fixed  in  the  geomagnetic  meridian 
(Vondrak  and  Baron,  1976).  Such  a scan  results  in  the  determination  of 
the  latitudinal  variation  of  measured  quantities.  The  elevation-scan 
data  used  in  this  study  were  slow  scans  from  the  northern  horizon  to 
the  southern  horizon,  covering  a range  in  geomagnetic  latitude  of  62° 
to  68°  in  the  E region  and  56°  to  73°  in  the  F region.  Each  scan 
required  approximately  15  minutes,  and  a 15-s  integration  time  was  used 
for  each  density  measurement. 

2.2  Data  Used  in  This  Study 

This  study  is  based  on  the  analysis  of  21  data  sets.  Each  data  set 
covers  an  interval  of  approximately  24  hours,  so  that  each  represents 
a continuous  measurement  at  all  local  times.  Two  complete  days  of 
elevation-scan  data  were  used.  Because  they  were  the  only  existing 
24-hour  runs  of  elevation-scan  data  at  the  time  the  study  began,  the 
description  of  auroral  morphology  in  Section  3 is  supplemented  by  a 
10-hour  run  on  19  February  1976.  In  addition,  19  runs  in  the  overhead 
mode  were  used.  Because  the  radar  sensitivity  and  data  quality  have 
been  continually  upgraded  in  recent  years,  data  sets  were  chosen  from  the 
time  interval  of  1974  to  1977,  The  data  are  listed  in  Table  2-1 
according  to  year  and  season.  A primary  criterion  used  in  data  selection 
was  to  have  data  representative  of  all  seasons  and  levels  of  magnetic 
activity.  As  can  be  seen  in  Table  2-1,  seven  data  sets  are  for  summer, 
five  dat.i  sets  are  for  winter,  and  nine  data  sets  are  for  equinox 
(including  spring  and  fall).  The  distribution,  with  respect  to  magnetic 
activity,  of  the  data  selected  for  this  study  is  shown  in  Figure  2-1. 

Kp  is  an  index  of  the  amount  of  global  magnetic  activity  during  a 3-hour 


Winter  Fall  Summer  Spring  I Season 


Table  2-1 


CHATANIKA  DATA  USED  IN  THIS  STUDY 


1974 

Date 

Time,  UT 

13  Feb 

0430  - 1930 

20  Mar 

0000  - 2400 

1975 

Date 

Time 

12  Mar 

0400  - 

15  Apr 

0130  - 

14  May 

0230  - 

17  Jul 

0820  - 

13  Aug 

0000  - 

17  Sep 

0000  - 

15  Oct 

0500  - 

17  Dec 

0140  - 

1976 


UT  Date 

2500 

17  Mar 

14  Apr 

1100 

12  May 

2400 

23  Jun 
14  Jul 

2400 

11  Aug 

2400 

2350 

2450 

20  Oct 
17  Nov 

2500 

22  Dec 

13  Nov* 
17  Dec* 


Time,  UT 

0030  - 2410 

0030  - 2350 

0000  - 2330 

0000  - 2400 
0000  - 2400 

0000  i 2330 

0010  - 2400 
0730  - 2410 

0020  - 2350 
0000  - 2400 
0000  - 2400 


* j 

Elevation-scan  data  i 

i 

1 
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interval.  For  comparison,  also  indicated  in  Figure  2-1  is  the  occurrence 
frequency  distribution  of  during  the  period  of  1932  to  1971  as  given 
by  Cage  and  Zawalick  (1972),  It  can  be  seen  that  the  data  selected  for 
this  study  have  been  chosen  to  represent  reasonably  well  the  distribution 
normally  occurring.  As  a result  we  would  expect  that  the  model  based  on 
these  data  would  have  the  minimum  statistical  error  for  conditions 
occurring  most  frequently. 


F«GURE  2-1  FREQUENCY  DISTRIBUTION  OF  Kp  FOR  THE  DATA  SET  USED 
IN  THIS  STUDY  AND  FOR  THE  YEARS  1932  TO  1971 
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3.  MORPHOLOGY  OF  THE  HIGH-LATITUDE  IONOSPHERE 

In  this  section  we  describe  the  morphology  of  the  high-latitude 
ionosphere  on  the  basis  of  Chatanika  incoherent-scatter  radar  data. 
Elevation-scan  (or  elscan)  data,  which  contain  electron-density-profile 
information  as  a function  of  both  latitude  and  local  time,  are  used  to 
identify  major  features  of  the  auroral  ionosphere  and  to  detect  any 
latitudinal  gradients  in  ionization.  Because  the  elscan  data  base  is 
so  limited,  we  supplement  it  with  overhead  radar  measurements  that  we 
used  to  examine  seasonal  and  diurnal  variations. 

The  characteristics  of  the  auroral  E region  are  clearly  evident  in 
the  data  presented  here.  Examination  of  the  data  indicated  that  the 
auroral  (particle  produced)  F layer  is  not  as  predictable  as  the  auroral 
E region.  With  the  data  base  available  for  this  study,  a substantial 
Fg  layer  was  found  to  be  present  at  night  only  Infrequently.  The  solar- 
controlled  F layer  dominated  during  the  day.  Consequently,  our  under- 
standing of  the  occurrence  of  a valley  region  in  the  dark  sector  is 
limited  by  the  unpredictable  nature  of  the  particle-produced  auroral 
Fg  layer. 

It  should  again  be  noted  that  the  electron  densities  presented 
here  have  not  been  corrected  for  either  the  temperature  ratio  or  the 
Debye-length  paramettr.  The  net  result  is  always  to  underestimate  the 
true  electron  density,  particularly  in  the  F region.  The  effects  of 
neglecting  the  altitude-dependent  correction  factor  are  discussed  in 
Appendix  A. 

3,1  Latitudinal  Variations  Evident  in  Elevation-Scan  Data 

Elscan  data  are  ideal  for  use  as  a basis  of  constructing  a 
morphological  description  of  the  auroral  ionosphere.  They  can  be  used 
to  construct  contour  maps  of  electron  density  in  coordinates  of  altitude 
and  latitude,  which  immediately  reveal  various  latitudinal  boundaries  in 


1 


23 


auroral  lonosphorlc  structure.  In  the  E region,  these  boundaries  include 
those  of  the  discrete  and  diffuse  aurora.  These  boundaries  can  be  related 
to  any  F-region  feature,  such  as  the  trough  region.  In  addition  to  their 
use  in  locating  latitudinal  boundaries,  the  contour  maps  may  be  used  to 
extract  any  latitudinal  gradients  that  may  exist  In  the  peak  electron 
density  and  Its  altitude.  By  collecting  elscan  data  continuously  over  a 
24-hour  period,  variations  in  the  above  characteristics  can  be  determined 
as  a fvinction  of  local  time  and  magnetic  activity. 

A description  of  high-latitude  morphology  from  these  data  has  two 
principal  limitations.  First,  an  elscan  with  the  Chatanika  radar  inter- 
sects the  E-region  ionosphere  from  62°  to  68°  geomagnetic  latitude,  or  a 
latitudinal  extent  of  6°.  Consequently,  we  can  characterize  the  aurora 
only  when  it  is  located  within  the  scanned  region.  Second,  the 
morphological  description  presented  in  this  section  is  based  on  only 
three  data  sets,  two  of  which  wer«  continuous  24-hour  runs.  Clearly, 
much  more  data  must  be  examined  under  a variety  of  conditions  before  all 
the  major  features  can  be  identified. 

The  elscan  data  used  to  construct  the  preliminary  model  are 
summarized  in  Table  3-1.  They  consist  of  two  continuous  24-hour  runs,  each 
beginning  and  ending  near  local  noon  (2200  UT) , and  a "partial"  run  that 
started  in  the  afternoon  and  ended  shortly  after  local  midnight  (1000  UT) . 
All  except  the  December  1976  data  set  were  taken  during  moderately  active 
geomagnetic  conditions  (i.e.,  the  index  averaged  a value  of  4 during 
the  ent j re  data-gathering  period) . 

T1 o local  geomagnetic  conditions  that  prevailed  during  the  three 
days  of  Interest  can  be  inferred  from  the  horizontal  perturbation  (H) 
component  of  the  geomagnetic  field  measured  at  College,  Alaska,  as  shown 
in  Figure  3-1.  The  magnitude  of  the  perturbations  can  be  estimated  from 
the  magnetograms  by  measuring  the  deviation  of  the  traces  from  the 
horizontal  line  that  represents  the  ambient  level  of  the  geomagnetic 
field. 
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ELEVATION  SCAN  DATA  SUMMARY 


Date 

(UT) 

Start  Time 
(UT) 

Stop  Time 
(UT) 

ZK 

P 

12-13  Nov  1976 

2230 

2230 

33 

16-17  Dec  1976 

2200 

2200* 

17+ 

19  Feb  1976 

0030 

1045 

31- 

* 

Not  in  elscan  mode  from  0950  to  1120  UT 


I 

I 

A qualitative  feeling  for  the  differing  degrees  of  magnetic 
■ - activity  during  the  three  days  can  be  seen  by  comparing  the  magnetograms 

in  Figure  3-1.  The  December  magnetogram  is  virtually  a straight  line, 
indicating  essentially  no  magnetic  activity  until  around  1500  UT  when  a 
200-Y  negative  bay  occurred.  In  contrast,  the  November  magnetogram  shows 
continuous  activity  throughout  the  24-hour  period  shown,  with  nearly 
100-Y  positive  bay  activity  from  2300  to  0600  UT,  and  a major  negative 
bay  (~  1300  Y)  starting  around  1030  UT.  The  February  magnetogram 
shows  conditions  that  are  somewhat  less  disturbed  than  those  shown  in 
the  November  magnetogram.  There  was  virtually  no  positive-bay  activity 
until  0500  UT.  The  largest  negative-bay  activity,  which  commenced  around 
120(1  UT,  was  weaker  than  the  major  bay  in  the  November  magnetogram  and 
occurred  after  1045  UT,  the  end  of  Chatanlka  radar  data  collection. 

In  the  following  subsections,  we  discuss  characteristics  of 
the  elscan  data  from  four  different  time  sectors:  day,  evening,  midnight, 
and  morning.  Each  time  sector,  as  defined  here,  can  be  thought  of  as  a 
time  period  (of  about  6 hours  duration)  centered  approximately  around 
, 1200,  1800,  0000,  and  0600  magnetic  local  times  (MLT). 

3.1.1  Day-Sector  Characteristics 

Day-sector  characteristics  of  electron  density  variations  with 
altitude  and  latitude  are  illustrated  in  Figure  3-2.  The  contour  plot 

L. 
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COLLEGE.  ALASKA 
12  NOVEMBER  1976 


13  NOVEMBER  1976 


16  DECEMBER  1976 


18  FEBRUARY  1976 


17  DECEMBER  1976 


19  FEBRUARY  1976 


UNIVERSAL  TIME 

I IGURE  3-1  COLLEGE  MAGNETOGRAMS  (horizontal  component  only)  SHOWING  MAGNETIC 
ACTIVITY  vs  TIME  FOR  THE  THREE  ELSCAN-RUN  DAYS.  The  unshaded 
portions  of  the  magnetograms  correspond  to  times  when  elevation-scan  data  were 
collected. 


FIGURE  3-2  CONTOURS  OF  EQUAL  ELECTRON  DENSITY  MEASURED  BY  THE  CHATANIKA 
RADAR  BETWEEN  2323  UT  AND  2336  UT  ON  13  NOVEMBER  1976  {1323  to 
1336  local  standard  time).  Contours  are  in  units  of  lO^/cm^. 


was  constructed  from  elscan  data  collected  around  magnetic  noon  on 
13  November  1976  when  the  solar  zenith  angle  over  Chatanlka  was  86° . 

The  geomagnetic  K Index  was  4+. 

P 

Under  quiet  geomagnetic  conditions,  the  day  sector  Is 
primarily  characterized  by  an  electron  density  distribution  that  Is 
controlled  by  solar  radiation.  As  Is  well  known,  and  as  Illustrated 
In  Figure  3-2,  under  these  conditions  the  Ionosphere  Is  vertically 
stratified  (except  around  sunset  and  sunrise)  with  a deterministic 
ele<;tron  density  profile  that  Is  closely  related  to  the  solar  zenith 
angle.  In  Section  3.2,  the  solar  control  of  the  electron  density  profile 
will  be  discussed  In  detail  using  overhead  data. 

When  there  Is  geomagnetic  activity,  an  E layer  produced  by 
particle  precipitation  often  develops.  The  extent  to  which  the  solar- 
produced  electron  density  profile  Is  modified  by  the  precipitation  source 
depends  on  the  prevailing  solar  zenith  angle.  In  the  example  shown 
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(Figure  3-2),  the  particle -produced  E layer  Is  clearly  evident.  The 
details  of  the  latitudinal  variation  are  shown  in  Figure  3-3,  in  which  we 
have  superimposed  vertical  profiles  extracted  at  40-km  Intervals  from  the 
contour  plot  in  Figure  3-2.  The  dense  E layer  present  to  the  north  of 
Chatanlka  is  seen  to  be  associated  with  a deep  valley  between  the  E and  F 
regions.  As  shown  later  with  the  overhead  data  at  these  solar  zenith 

t 

angles,  a deep  valley  does  not  develop  when  solar  illumination  is  the 
only  source  of  ionization.  Therefore,  the  E layer  that  varies  with 
latltixle  and  the  corresponding  valley  between  the  E and  F regions  are 

both  controlled  by  particle  precipitation.  ] 


ELECTRON  DENSITY  — cm"® 


FIGURE  3-3  ELECTRON  DENSITY  PROFILES  EXTRACTED  AT  40-km 
INTERVALS  FROM  THE  CONTOUR  PLOT  IN  FIGURE  3-2 
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The  depth  of  the  valley  is  dependent  on  the  pcik  E-layer 
electron  density  produced  by  particle  pre<  Ipltation,  the  altitudes  of 
the  E and  F layers,  and  the  F^^-layer  electron  density  which  is  under 
the  control  of  the  solar  zenith  angle.  During  summer,  when  the  solar 
zenith  angle  is  much  smaller  for  a given  local  time,  the  valley  is  largely 
"filled  in"  by  the  presence  of  significant  F^^-layer  electron  density.  In 
this  til  e sector,  the  altitude  of  the  peak  E-layer  electron  density  is 
typically  100  km  or  lower,  which  allows  for  the  presence  of  a deeper 
valley  than  in  other  time  sectors.  The  relatively  low  peak  altitude  of 
the  E layer  suggests  that  the  morphology  of  the  related  precipitation 
may  be  evident  in  the  morphology  of  riometer  absorption. 

The  presence  of  a day-sector  particle-produced  E layer  near 
Chataiiika  (which  is  located  at  geomagnetic  latitude)  is  consistent 
with  1 he  findings  of  other  workers.  Plddlngton  (1965)  and  Hartz  and 
Brice  (1967)  organized  a large  body  of  data  related  to  particle  preci- 
pitation and  concluded  that  there  are  two  zones  of  auroral  particle 
precipitation.  One  essentially  coincides  with  the  auroral  oval  that 
contains  the  visual  aurora  and  i.s  far  poleward  of  Chatanlka  during  th<‘ 
day.  The  other  is  a ring  of  precipitation  centered  at  65°  geomagnetic 
latitude.  Independent  of  local  time.  Hartz  and  Brice  '1967)  found  that  ■! 

the  occurrence  maximum  for  the  second  zone  (their  "drl/.zle"  zone)  was  i 

around  0800  MLT.  The  "drizzle  zone"  precipitation  has  also  been  reported  j 

from  a 11-sky -camera  observations  (Davis , 1962)  and  ISIS-2  scanning  ? 

photometer  observations  (Lui  et  al.,  1973).  It  is  responsible  for  one  of  ; 

the  characteristic  diurnal  maxima  in  auroral  absorption,  as  discussed  in  ' 

Section  6. 

The  only  occasion  when  large  latitudinal  gradients  are  observed 
in  the  dayside  F region  is  at  sunset  and  sunrise.  Data  taken  during  an 
elevation  scan  at  local  sunset  at  Chatanlka  are  shown  in  Figure  3-4.  At 
this  time  the  sun  has  set  in  the  north  but  is  still  illuminating  the 
F region  at  more  southerly  latitudes.  The  effects  of  sunset  are  very 
apparent.  In  the  north  the  maximvim  electron  density  has  decreased  and 
the  Isodensity  contours  are  tilted  upward  as  the  lower  F-reglon  erodes 
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19  FEB  1976 
0221-0236  UT 

FIGURE  3-4  CONTOUR  MAP  OF  EQUAL  ELECTRON  DENSITY  MADE  FROM  CHATANIKA 
ELEVATION-SCAN  DATA  TAKEN  AROUND  LOCAL  SUNSET.  A sunlit  F region 
it  seen  to  the  south  of  Chatanika  white  the  post-sunset  decay  of  the  F region  is 
evident  to  the  north. 


due  to  faster  recombination  at  lower  altitudes.  In  the  south  the  sun  Is 
still  illuminating  the  ionosphere  and  the  F-reglon  Is  horizontally 
stratified  similar  to  the  sunlit  ionosphere  shown  in  Figure  3-2. 

3.1.2  Evening-Sector  Characteristics 

A contour  plot  of  electron  density  representative  of  the 
evening  sector  is  shown  in  Figure  3-5.  The  data  set  was  taken  around 
1930  ML'  on  19  February  1976  when  the  solar  zenith  angle  at  Chatanika 
was  11'/  . Since  the  solar-controlled  E layer  is  negligible  for  solar 
zenith  angles  greater  than  90°  (see  Watt , 1973;  and  next  section),  the 
electron  density  distribution  shown  in  Figure  3-5  is  due  solely  to 
particle  precipitation.  The  geomagnetic  K index  was  4. 
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The  contour  plot  In  Figure  3-5  can  be  chiiracterlzed  by  a 

region  of  relatively  uniform  particle  precipitation  centered  on  Chatanika 

and  extending  ±200  km  in  latitude,  and  by  an  Isolated  but  more  intense 

particle  precipitation  region  centered  270  km  north  of  Chatanika.  The 

first  rescion  is  associated  with  the  diffuse  visual  aurora  and  is 

characterized  by  a relatively  uniform  peak  electron  density  of  approxi- 
5 3 

mately  2 x 10  el/cm  and  located  at  an  altitude  of  120  km.  The  second 

region,  which  is  associated  with  the  discrete  visual  aurora,  is 

6 3 

characterized  by  a peak  electron  density  approaching  10  el/cm  with  its 
altitude  located  closer  to  105  km.  The  more  structured  contours,  which 
are  associated  with  t)ie  discrete  visual  aurora,  are  generally  located 
near  the  poleward  boundary  of  the  region  of  diffuse  aurora. 

The  structure  in  the  diffuse  aurora  (e.g.,  that  seen  over 
Chatanika  in  Figure  3-5)  is  probably  related  to  substorm  activity. 

During  less  disturbed  periods,  the  topside  of  the  E layer  is  often  very 
uniform,  similar  in  character  to  that  shown  in  Figure  3-2. 

The  F-reglon  electron  densities  in  Figure  3-5  are  much 
smaller  than  those  in  the  auroral  E layer.  This  feature  appears  to  be 
representative  of  the  evening  sector,  as  well  as  the  remaining  dark 
sector  (see  following  sections).  The  decrease  in  electron  density 
with  altitude  is  usually  monotonic  so  that  no  valley  is  present.  The 
only  conspicuous  evening  F-region  feature  is  an  enhancement  occasionally 
found  at  the  equatorial  edge  of  the  diffuse  aurora.  Such  an  enhancement 
is  visible  in  Figure  3-5  at  a distance  of  150  km  south  of  Chatanika. 

This  enhancement,  referred  to  as  the  "boundary  feature,"  is  generally 
transient  in  nature,  persisting  for  typically  30  minutes. 

3.1.3  Midnight-Sector  Characteristics 

The  characteristics  of  the  latitudinal  distribution  of  electron 
density  in  the  midnight  sector  are  not  unlike  those  seen  in  the  evening 
sector.  The  latitudinal  distribution  is,  however,  much  more  complex  and 
dynamic  In  nature,  a behavior  directly  related  to  substorm  effects.  In 
the  ever  Lng  sector,  virtually  the  only  substorm  effects  are  those  that 
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FIGURE  3-5  ISODENSITY  CONTOUR  MAP  OF  THE  EVENING-SECTOR  AURORAL 
IONOSPHERE.  The  auroral  E layer  over  Chatanika  corresponds  to  the 
diffuse  aurora,  and  the  localized  structure  at  270  km  north  of  Chatanika 
corresponds  to  the  discrete  aurora.  The  F-region  boundary  feature  is 
seen  near  the  equatorward  boundary  of  the  diffuse  aurora. 

propagati  from  the  midnight  sector,  behind  the  westward-traveling  surge. 
These  effects  are  usually  contained  along  the  poleward  boundary  of  the 
diffuse  aurora.  In  the  midnight  sector,  the  diffuse  aurora  appears  to 
bo  directly  related  to  the  substorm  process  that  originates  there. 

A contour  plot  illustrating  the  E-layer  electron  density  distri- 
bution with  latitude  in  the  midnight  sector  is  shown  in  Figure  3-6.  The 
data  were  collected  around  2230  MLT  on  13  November  1976,  when  the  solar 
zenith  angle  was  132°.  Referring  to  Figure  3-1,  we  see  that  the  magnetic 
H component  had  Just  returned  to  the  quiescent  level  from  positive-bay 
activity.  While  these  data  should  not  be  interpreted  as  data  representa- 
tive of  quiet  times,  they  may  be  interpreted  as  representative  of  non- 
substorm conditions. 
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CONDITIONS  ON  13  NOVEMBER  1976 


The  electron  density  contours  In  Figure  3-6  Indicate  the 
presence  of  a relatively  uniform  diffuse  aurora,  extending  from  1°  south 
to  3°  north  of  Chatanlka . Peak  E-layer  electron  density  Is  comparable 
to  that  found  In  the  evening  sector  (see  Figure  3-5).  The  only  noticeable 
difference  Is  a slight  lowering  of  the  altitude  of  peak  electron  density 
to  approximately  115  km.  As  In  the  evening  sector,  the  electron  densities 
In  the  F region  are  much  smaller  than  those  In  the  E layer.  We  show 
later  In  this  section  that  the  decrease  Is  essentially  monctonlc  with 
altitude,  with  no  evidence  of  a significant  valley. 

Referring  to  Flgui-e  3-1,  we  note  the  onset  of  a substorm  around 
1030  ITT  on  13  November  1976  with  continuous  negative-bay  activity  there- 
after. I'o  compare  the  behavior  of  the  diffuse  aurora  in  the  midnight 
sector  under  substorm  and  non-substorm  conditions,  we  have  shown  in 
Figure  3-7  a contour  plot  constructed  from  data  taken  around  2330  MLT 
(1100  UT)  on  the  same  day.  The  diffuse  aurora  In  Figure  3-7  is  seen  to 
extend  over  a much  larger  latitudinal  extent  than  that  scanned  by  the 
radar  (€•  ).  The  peak  electron  density  Is  seen  to  be  much  greater  than 
that  in  I'igure  3-6  at  an  altitude  closer  to  100  km  than  to  115  km.  The 
base  of  the  Ionosphere  is  found  in  Figure  3-7  to  be  at  nearly  80  km  as 
compared  to  100  km  in  Figure  3-6. 

While  the  electron  densities  in  the  E layer  are  found  to  be 

consider  ibly  enhanced  during  substorm  activity,  there  Is  no  apparent 

change  i i the  F region.  The  F-region  electron  densities  in  Figure  3-7 

do  not  d ffer  significantly  from  those  in  Figure  3-6,  remaining  below 
5 3 5 3 

10  el/cm  . While  the  10  -el/cm  contour  Is  seen  to  extend  up  to  nearly 
230  km  altitude.  It  may  be  thought  of  as  representing  the  topside  of  the 
E layer.  In  any  case,  a valley  is  not  a prominent  nor  persistent  feature 
in  the  midnight  sector,  under  either  substorm  or  non-substorm  conditions. 

3.1  4 Morning-Sector  Characteristics 

The  diffuse  aurora  is  also  found  in  the  morning  sector.  However, 
according  to  Akasof u (1977),  the  diffuse  aurora  in  this  time  sector  is 
characterized  by  nonuniform  luminosity.  In  contrast  to  the  relatively 
structureless  nature  of  the  evening  diffuse  aurora. 
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A contour  plot  of  the  diffuse  aurora  that  occurred  around 

0320  MLT  on  13  November  1976  Is  shown  In  Figure  3-8.  This  plot,  from  a 

period  of  recovery  from  negative-bay  activity  (see  Figure  3-lX  is  seen 

to  be  similar  to  that  found  in  the  midnight  sector.  The  peak  electron 

5 3 

density  is  much  greater  than  2 x 10  el/cm  and  its  altitude  is  around 
100  km.  The  latitudinal  extent  of  the  precipitation  region  appears  to 
be  comparable  to  or  greater  than  the  6°  latitudinal  sector  scanned  by  the 
radar.  As  before,  the  electron  density  decreases  monotonies lly  with 
altitude  above  the  E layer,  with  no  valley  between  the  E and  F regions. 
This  appears  to  be  characteristic  of  the  dark  sector. 

A second  contour  plot  constructed  from  data  taken  later  that 

day  (0555  MLT),  during  a period  of  relatively  weak  magnetic  activity,  is 

shown  in  Figure  3-9.  At  this  time,  the  solar  zenith  angle  over  Chatanika 

was  95° , The  diffuse  aurora  is  still  widespread,  although  somewhat  patchy 

over  Chatanika  The  peak  electron  density  has  decreased  to  less  than 
5 3 

2 X 10  el/cm  but  its  altitude  has  decreased  to  around  95  km.  A close 
examination  of  Figure  3-9  also  reveals  that  the  altitude  of  the  E-layer 
peak  increases  somewhat  with  decreasing  latitude. 

3,1,5  Summary  of  Diurnal  Variations  Deduced  from  Elscan  Data 

In  order  to  Illustrate  the  characteristic  profile  features  more 
clearly,  we  have  taken  the  12-13  November  1976  data  set  and  replotted 
selected  elscans  as  sets  of  superimposed  electron  density  profiles.  The 
results  are  shown  in  Figure  3-10. 

The  superimposed  electron  density  profiles  shown  for  a given 
time  were  computed  from  elscan  data  at  40-km  intervals  over  ±240  km  extent 
centered  over  Chatanika.  Arrows  from  each  set  of  profiles  point  to 
corresponding  times  on  a plot  using  geomagnetic  coordinates  of  latitude 
vs  local  time.  The  "clock  dial"  plot  in  the  center  of  Figure  3-10  shows 
(1)  the  latitude  of  Chatanika  by  a dashed  circle  at  65°,  (2)  the 

5 3 

latitudinal  extent  of  E-layer  electron  densities  greater  than  10  el/cm 
for  each  elscan  by  radial  line  segments,  and  (3)  the  location  of  the 
auroral  oval  for  a magnetic  Q index  of  5 by  the  shaded  region.  Sunrise- 
sunset  lines  for  altitudes  of  100  and  300  km  also  are  shown. 
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PERIOD  FROM  NEGATIVE-BAY  ACTIVITY  ON  13  NOVEMBER  1976 


ELECTRON  DENSITY  — 10  cm 


RELATIVELY  WEAK  MAGNETIC  ACTIVITY  ON  13  NOVEMBER  1976 


The  day  sector  is  characterized  by  the  profiles  shown  for 
times  1110,  1302,  and  1501  MLT.  We  note  the  presence  of  a stable 
stratified  F region  at  1110  MLT  when  the  solar  zenith  angle  was  near 
its  minimum  diurnal  value  of  83° . The  progressive  decay  of  the  F layer 
is  seen  in  the  profiles  for  1302  and  1501  MLT.  The  rapid  decay  of  the 
solar-controlled  F^^  layer  is  an  Important  feature  of  the  day  sector. 

Auroral  precipitation  during  magnetically  active  conditions,  as  shown  in 
this  data  set,  results  in  the  development  of  a pronounced  valley.  The 
variability  in  the  E-layer  electron  density  profiles  is  due  primarily  to 
latitudinal  gradients  and  structure,  although  rapid  temporal  variations 
can  coifuse  the  observed  spatial  structure,  especially  during  periods  of 
substorm  activity. 

The  evening  sector  is  characterized  by  profiles  shown  for  times 
1701,  1906,  and  2058  MLT.  Transition  from  day  sector  to  evening  sector 
characteristics  is  seen  in  the  profiles  for  1701  MLT.  Because  the  F- 
layer  electron  densities  have  decayed  to  values  difficult  to  measure 
accurately  in  the  elscan  mode,  the  apparent  structure  at  those  altitudes 
shoult'  not  be  Interpre  ed  as  actual  electron  density  variations.  The 
low-altitude  E layer  s<‘en  in  one  of  the  profiles  was  found  at  the  lowest 
latitude.  The  remaining  profiles  show  a progressive  buildup  in  peak 
electron  density  with  Increasing  latitude.  Except  for  the  profile  with 
a low  E layer,  which  is  similar  to  those  seen  in  the  day  sector,  the 
remaining  profiles  can  be  shown  to  have  properties  corresponding  to  the 
diffuse  aurora.  All  have  nearly  Identical  peak  altitudes  and  slopes  of 
the  topside  of  the  E layer. 

Examining  the  profiles  from  the  1906  and  2058  MLT,  we  find 
that  the  profile  shapes  are  generally  similar  except  for  a progressive 
lowering  of  the  peak  altitude  with  increasing  local  time,  and  the 
apparent  transitory  appearance  of  a valley  region  at  1906  MLT.  Noting 
that  the  gradient  on  the  topside  of  the  E layer  remains  virtually  constant, 
it  appears  that  the  existence  of  a valley  region  is  directly  dependent  on 
the  development  of  a significant  F^  layer.  If  the  electron  density 
profiles  above  the  E layer  observed  around  1906  MLT  are  due  to  particle 
preclpil itlon,  the  Incoming  particles  may  have  a non-monotonlc  energy 
spectrum 
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The  midnight -sector  characteristics,  illustrated  by  the 
profiles  for  0002  MLT,  do  not  appear  to  differ  much  from  the  profiles 
from  the  evening  sector.  There  is  a tendency  for  the  topside  slope  to 
flatten  around  200  km  altitude,  suggesting  some  enhancement  in  the  F 
region.  However,  examination  of  the  other  elscans  did  not  reveal  a 
persistent  F layer.  The  two  profiles  with  the  greatest  peak  electron 
densities  occurred  at  the  two  extremes  (north  and  south)  of  the 
latitudinal  scan.  Both  appear  to  be  associated  with  discrete  visual 
aurora . 

The  morning-sector  characteristics  are  illustrated  by  profiles 
from  0207,  0416,  and  0629  MLT.  Perhaps  the  most  apparent  feature  in  the 
profiles  from  this  sector  is  the  relatively  large  peak  electron  density 
of  the  E layer  with  steeper  topside  slope  than  in  other  time  sectors.  As 
mentioned  previously,  the  E-layer  peak  altitude  is  much  lower  than  that 
in  the  evening  and  midnight  sectors.  Another  characteristic  of  the 
morning  E layer  is  the  lowering  of  the  peak  altitude  with  decreasing 
latitude,  as  seen  in  the  profiles  from  0207  MLT.  Again,  there  is  an 
absence  of  a persistent  F layer,  and  therefore  no  valley  region.  The 
appearance  of  a valley  is  seen  in  the  0629  MLT  profiles  with  the  develop- 
ment of  a sunlit  Fg  layer. 

3.2  Diurnal  Variations  Evident  In  Overhead  Data 

The  elevation-scan  data  shown  in  the  preceding  section  are  all  from 
the  winter  season.  In  this  section  we  will  assess  the  deterministic 
characteristics  of  the  solar-controlled  electron  density  profile  by 
showing  overhead  data  from  other  seasons. 

We  have  selected  three  24-hour  runs  as  examples.  These  data  sets 
are  summ  irized  in  Table  3-2.  The  data  sets  were  selected  primarily 
because  of  the  relatively  quiet  magnetic  conditions  that  prevailed  over 
each  24-hour  period.  Quiet  conditions  were  desired  so  that  solar  control 
of  the  ionospheric  layers  could  be  determined  with  a minimum  amount  of 
contamination  by  particle  precipitation. 
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Table  3-2 


OVERHEAD  DATA  SUMMARY 


Date 

(UT) 

Start  Time 
(UT) 

Stop  Time 
(UT) 

T.K 

P 

23  June  1976 

0100 

2400 

8+ 

17  December  1975 

0200 

2400 

16- 

17  March  1976 

0100 

2400 

27- 

Figure  3-11  shows  the  magnetograms  of  the  horizontal  component  of 
the  geomagnetic  field  as  measured  at  College  during  the  three  days  of 
Interest.  We  note  that  the  magnetic  perturbations  on  summer  and  winter 
days  are  comparable,  and  very  similar  to  that  presented  in  the  previous 
section  for  the  17  December  1976  elscan  data  set.  The  magnetogram  for 
the  equinox  example  (17  March  1976)  is  seen  to  be  more  active  than  the 
other  two,  and  similar  to  the  magnetic  conditions  that  existed  during 
the  19  February  1976  elscan  data  set  (see  Figure  3-1) . 

The  electron  density  profiles  from  the  three  24-hour  runs  are 

presented  in  Figure  3-12.  Each  row  contains  approximately  one-half  of 

the  entire  data  set  for  that  day.  For  each  hour,  six  individual  profiles 

are  superimposed.  Each  represents  a five-minute  average  taken  within  15 

minutes  of  the  hour.  The  solar  zenith  angle  and  K index  for  each  hour 

P 

are  also  shown.  In  collecting  the  data,  the  radar  was  operated  in  a 
three -posit  ion  dlscrete-azitquth  scan  with  a fixed  elevation  angle  of 
76.5° . Therefore,  some  of  the  variability  in  the  slopes  of  the  super- 
imposed profiles  is  due  to  horizontal  gradients  and  structure  in  the 
ionosphere . 

let  us  consider  first  the  behavior  of  the  summer  ionosphere  above 
Chatanika  over  a 24-hour  period.  On  23  June  1976,  the  solar  zenith 
angle  was  always  less  than  91° . Therefore,  the  F^  region  was  always 
sunlit.  The  profiles,  starting  at  0100  UT  (1500  LT) , are  seen  to  have 
a very  regular  shape  which  persists  until  0500  UT  when  the  effects  of 
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FIGURE  3-12  ELECTRON  DENSITY  PROFILES  FROM  ±15  MINUTES  ABOUT  THE  HOUR  FOR  THREE  OVERHEAD 
24-HOUR  RUNS.  Each  set  of  profiles  shown  on  the  hour  consists  of  approximately  six  individual  profiles, 
each  averaged  over  five  minutes. 


sunset  in  the  E and  region  began.  The  decay  of  the  two  ionospheric 
layers  is  complete  by  1000  IIT.  In  the  local  midnight  and  early-morning 
sectors,  we  see  the  only  evidence  of  a particle-produced  E layer.  A 
valley  appears  only  around  1100  UT  when  the  solar  zenith  angle  is 
greatest  (and  hence,  corresponds  to  a time  when  the  electron  density 
is  minimum).  From  1200  UT  on,  we  see  the  sunrise  effects  in  the  E and 
Fj^  regions.  The  period  between  1200  and  1400  UT  contains  nighttime 
profiles  without  an  auroral  E layer  slightly  south  of  Chatanlka , and 
profiles  with  a superimposed  auroral  E layer  slightly  north  of  Chatanlka. 
This  indicates  that  the  equatorward  boundary  of  the  diffuse  aurora  was 
centered  over  Chatanlka.  By  1500  UT,  the  diffuse  aurora  had  moved  pole- 
ward  and  we  see  electron  density  profiles  controlled  by  solar  illumination. 
The  profiles  from  1600  UT  to  2400  UT  are  essentially  replicas  of  the  first 
profile  shown  at  0100  UT. 

We  can  now  compare  the  diurnal  variations  of  a quiet  sununer  ionosphere 
to  those  of  a quiet  winter  ionosphere,  shown  in  the  second  row  of 
Figure  3-12.  The  winter  profiles  around  0200  IT  display  a well  behaved 
Fg  layer  and  a low-altitude  auroral  E layer.  If  we  scan  the  summer  and 
winter  profiles  (first  and  second  rows  in  Figure  3-12)  for  similar  F- 
layer  shapes,  we  find  that  the  solar  zenith  angles  are  similar.  The 
major  difference  is  a slightly  larger  peak  Fg-layer  electron  density 
during  the  winter.  The  winter  auroral  E layer  is  similar  to  that 
seen  at  the  same  local  time  and  season  under  more  disturbed  magnetic 
conditions  (see  Figure  3-10) . The  decay  of  the  layer  with  increasing 
solar  zenith  angle  is  seen  to  occur  between  0300  and  0500  UT.  There  is 
evidence  in  the  profiles  around  0700  UT  that  a measurable  Fg  layer,  one 
that  is  not  solar-controlled,  can  occasionally  appear.  But  as  we  found 
in  the  previous  section  (and  through  examination  of  all  the  overhead 
data),  it  is  not  a persistent  nor  reliable  feature  of  the  nighttime 
ionosphere. 

At  0800  UT,  a weak  auroral  E layer  is  seen  at  an  altitude  that  is 
measurably  higher  than  that  of  the  E layer  seen  during  the  previous  hour. 
The  higher-altitude  E layer  seen  between  0800  and  1300  UT  is  the  evening 
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and  midnight  diffuse  aurora.  The  characteristic  lowering  of  the  peak 
altitude  of  the  auroral  E layer  in  the  morning  sector  is  seen  between 
1400  and  1600  MT. 

The  sunrise  effects  can  be  seen  between  1800  and  2000  UT.  Because 
of  the  large  zenith  angles,  the  E and  layers  never  develop  fully. 
Instead,  they  form  a ledge  of  constant  electron  density  up  to  the  F 

A 

layer.  The  Fg  layer  profiles  between  2000  and  2400  UT  are  similar  to 
the  profile  around  0900  UT  in  the  first  row,  except  for  a slightly 
steeper  topside  slope.  Variability  in  the  topside  slope  may  be  due  to 
temperature  effects  in  the  incoherent -scatter  data,  which  have  not  been 
corrected  for  in  these  results. 

Finally,  we  examine  the  profiles  shown  in  the  third  row  in  Figure 
3-12,  taken  during  an  equin.ictial  period,  under  more  disturbed  magnetic 
conditions.  The  sunset  behavior  seen  in  the  summer  ionosphere  is  again 
evident  in  this  data  set,  except  that  it  occi'.rs  at  an  earlier  local  time 
(but  at  similar  solar  zenith  angles).  The  low-altitude  auroral  E layer 
is  also  evident  from  0200  to  0600  UT,  when  the  evening  diffuse  aurora 
entered  the  Chatanika  radar  beam.  In  this  data  set,  there  is  very  little 
evidence  of  even  a transient  layer  in  the  dark  sector.  The  dark 
sector  is  dominated  by  the  diffuse  aurora  (l.e.,  the  auroral  E layer). 

The  sunrise  effect  in  the  E region  is  masked  by  the  presence  of  the 
particle -produced  E layer.  However,  the  stable  solar-controlled 
ionosphere  is  evident  after  1900  UT. 

From  the  above  discussion,  it  is  clear  that  the  quiet  high-latitude 

electron  density  profile  is  controlled  by  the  solar  zenith  angle.  This 

control  is  most  clearly  seen  during  sunset  and  sunrise.  Electron  density 

profiles  that  Illustrate  the  decay  of  the  E and  F^  layers,  from  two  summer 

days,  are  shown  in  Figure  3-13.  For  solar  zenith  angles  less  than,  say, 

70° , the  electron  density  is  characterized  by  a nearly  flat  profile  with 

evidence  of  an  E layer  at  105  km,  an  F layer  at  150  km,  and  an  F_  layer 

1 » 

around  250  km.  Peak  electron  densities  for  all  three  layers  are  of  the 
5 3 

order  of  10  el/cm  . (We  should  note  again  that  these  profiles  have  not 
been  corrected  for  the  electron-to-ion  temperature  ratio,  T^.)  The  decay 
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of  the  E and  layers  during  sunset  is  seen  to  occur  as  the  solar  zenith 
angle  Increase  from  75°  to  96° . Such  a systematic  decrease  In  the  density 
at  120  km  has  been  shown  by  Watt  (1973).  The  presence  of  a small  E layer 
greater  than  the  layer  may  simply  be  due  to  weak  particle  precipitation 
Finally,  the  topside  scale  height  ranges  between  123  km  and  196  km,  If 
we  use  the  profiles  taken  on  13  May  1976. 
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4.  CHATANIKA  MODEL  OF  THE  HIGH-LATITUDE  IONOSPHERE 


In  this  section  we  describe  the  development  of  the  ionospheric 
model  derived  from  the  Chatanika  data.  This  model  is  basically  the 
RADC-F>0‘  AR  model  with  an  improved  specification  of  the  auroral  E layer 
and  the  valley  region  above  it. 

Like  the  RADC-POLAR  model,  our  modified  version  has  an  empirical 
basis.  Each  of  the  Chatanika  profiles  used  in  this  study  was  fitted 
parametrically.  The  derived  parameters  were  then  statistically  related 
to  variations  in  solar  zenith  angle  and  geomagnetic  activity. 

The  first  part  of  this  section  describes  the  parametric  fitting. 

Next  we  discuss  the  correlation  of  the  parameters  with  solar  and  geo- 
magnetic variations  and  the  types  of  profiles  constructed  in  the  model 
from  the  parameters.  This  is  followed  by  a discussion  of  the  latitude 
and  local-time  variation  in  the  auroral-E  layer  specified  in  the  model. 
Finally  we  compare  old  and  new  profile  predictions  with  Chatanika  data 
on  a specific  day. 

4 . 1 Extraction  of  Profile  Parameters  from  Chatanika  Data 

Ihe  ionosphere  generator  in  the  RADC-POLAR  model  employs  eight 
parameters  to  generate  electron  density  profiles  throughout  the  E,  F^, 
and  Fg  regions  at  a particular  point  in  space  and  time.  These  parameters 
are:  height  of  maximum  electron  density  of  the  E,  F , and  F„  layers; 
critical  frequency  of  the  three  layers;  and  semi-thickness  of  the  E and 
F^  layers. 

In  the  development  of  the  RADC-POLAR  model,  extensive  data  sets 
were  used  in  determining  the  morphology  of  the  eight  parameters.  In 
particular,  the  F-layer  specification  in  that  model  is  more  accurate 
than  could  be  derived  from  the  Chatanika  data  base  used  in  this  study. 
Therefore,  we  retained  the  parameters  of  the  Fg  layer.  Our  aim  in 
modifying  the  model  was  to  preserve  its  major  predictions  and  morphology 
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while  introducing  a more  realistic  profile  shape  and  Improved  auroral 
E parameters. 


r 


We  have  added  four  new  parameters  to  define  potential  valleys  in  the 
region  between  the  E and  F layers.  If  a single  minimum  Is  present,  it  is 
referred  to  as  valley  one.  If  an  layer  is  present,  the  minimum 
between  the  E and  F^  layers  is  referred  to  as  valley  one,  and  the  minimum 
between  the  F and  F layers  is  referred  to  as  valley  two.  The  four  new 

X M 

valley  parameters  are: 

(1)  foVl/fo^2.  the  ratio  of  the  plasma  frequency  of  valley  one 
to  the  critical  frequency  of  the  F2  layer. 

(2)  ^o'^2'^^o^2’  ratio  of  the  plasma  frequency  of  valley  two 

to  the  critical  frequency  of  the  layer. 

(3)  (HVj^  - HEl/CHFg  - HE),  the  height  above  the  E layer  of 
valley  one,  measured  as  a fraction  of  the  separation 
between  the  F^  layer  and  the  E layer. 

(4)  (HVg  - HEl/CHF^  - HE),  a similar  height  specification 
for  valley  two. 

It  should  be  noted  that  we  have  loosely  used  the  term  "f^^  layer"  to 
represent  any  local  maximum  between  the  E and  F^  layers;  this  layer 
is  not  necessarily  the  conventional  solar-produced  F^^  layer.  The  four 
valley  parameters  are  used  in  conjunction  with  the  original  eight 
parameters  to  produce  realistic  profiles  that  more  nearly  resemble  those 
found  in  the  Chatanika  data. 

Each  Chatanika  profile  consists  of  electron  density  measurements 
at  5-km  altitude  intervals.  The  profiles  from  the  overhead  data  were 
computed  from  integration  times  equal  to  either  the  dwell  time  at  each 
position  or  the  azimuthal  scan  time  (typically  5 or  6 minutes).  For  the 
elevation-scan  data  we  extracted  from  each  contour  the  vertical  profile 
overhead  at  Chatanika.  From  each  profile,  12  parameters  (8  original  and 
4 valley)  were  determined.  Since  the  number  of  profiles  was  quite  large, 
a method  for  finding  the  persmeters  was  needed  that  would  be  fast, 
efficient  and  consistent.  A computer  program  was  developed  that  found 
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the  12  parameters  automatically.  Improvements  were  made  on  this  program 
as  the  data  processing  was  done,  so  that  all  of  the  data  were  not 
processed  Identically.  However,  before  the  final  statistics  were  derived 
from  the  parameter  values,  all  of  tlie  data  were  checked  and  poorly  fitted 
values  were  eliminated.  The  final  data  base  was  a consistent  set  of 
parameters  that  were  then  used  to  determine  the  morphology  of  the  fitting 
parameters , 

4 . 2 Variation  of  Profile  Parameters  with  Solar  Zenith  Angle  and 

Geomagnetic  Activity 

The  morphology  of  the  parameters  has  been  investigated  in  terms  of 
magnetic  activity  and  solar  . Illumination,  the  phenomena  that  produce  the 
largest  effects  in  the  auroral  ionosphere.  The  two  categories  represent 
two  different  physical  phenomena:  (1)  the  changes  in  electron  density 
associated  with  magnetic  activity  general’ly  caused  by  auroral  particle 
precipitation,  and  (2)  the  changes  related  to  solar  photoionization. 


As  an  index  for  sorting  the  solar  effects  we  chose  solar  zenith 
angle,  which  allowed  us  to  combine  both  diurnal  and  seasonal  changes 
into  one  index.  The  data  were  grouped  into  five  categories  of  solar 
zenith  angle:  x < 60°  , 60°  < x < 80°  , 80°  < x < 90°  , 90°  < x <"  100°  , 
and  X ^ 100°  . 


As  an  index  for  magnetic  effects  it  was  desirable  to  use  K , since 

P 

that  parameter  is  readily  available  and  is  in  general  use  in  predictive 

programs.  is  a planetary  geomagnetic  activity  index  and  is  derived 

from  12  midlatitude  magnetic  observatories.  We  used  K as  a first  try 

P 

at  predicting  the  changes  in  electron  density  as  a function  of  magnetic 
activity.  However,  we  suspected  that  the  correlation  might  be  improved 
by  using  a hlgh-latltude  magnetic  index  rather  than  a planetary  one.  We 
chose  the  local  magnetic  index  at  College,  Alaska,  located  approximately 
30  km  from  Chatanika. 


Scatter  plots  of  parameter  values  as  a function  of  K and  as  a 

P 

function  of  K were  prepared.  Figures  4-1  and  4-2  show  such  scatter  plots 
for  one  of  the  parameters — maximum  electron  density  of  the  auroral  E 
layer — for  the  two  largest  solar-zenith-angle  categories.  This  parameter 
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FIGURE  4-1  CORRELATION  OF  THE  MAXIMUM  ELECTRON  DENSITY  OF  THE  AURORAL 
E LAYER  WITH  VARYING  GEOMAGNETIC  ACTIVITY  FOR  SOLAR  ZENITH 
ANGLES  OF  90-100°.  Th«  geomagnetic  indices  used  are  Kp  (upper  piot)  and 
K (iower  plot). 
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FIGURE  4-2  CORRELATION  OF  THE  MAXIMUM  ELECTRON  DENSITY  OF  THE 
AURORAL  E LAYER  WITH  VARYING  GEOMAGNETIC  ACTIVITY 
FOR  SOLAR  ZENITH  ANGLES  > 100°.  The  geomagnetic  indices 
used  are  Kp  (upper  plot)  and  K (lower  plot). 


is  the  residual  after  solar  E has  been  subtracted  from  the  data.  Note 
that  the  curves  appear  different  partly  because  finer  quantization 


Increments  are  available  for  the  K index.  However,  there  is  more  visual 

P 

evidence  of  a trend  with  the  K index  than  there  is  with  the  K index. 

P 

For  each  of  the  parameters  analyzed  in  the  statistical  program, 
linear  least-squares  fits  were  made  to  the  data  in  each  category  of 
solar  zenith  angle.  The  correlation  coefficients,  R,  and  number  of 
points  used,  N,  are  shown  In  Tables  4-1  and  4-2.  Also  in  the  tables  we 
have  listed  the  standard  error  of  the  estimate,  a,  which  is  one  standard 
deviation  from  the  estimate.  The  confidence  assumed  in  the  R and  a 
values  depends  on  the  number  of  points.  Correlations  computed  from  small 
samples  that  are  not  significant  at  the  95%  confidence  level  have  been 
marked  with  an  asterisk. 

The  six  parameters  listed  in  Tables  4-1  and  4-2  are  the  two 
E-layer  parameters  and  four  valley  parameters  that  have  been  incorporated 
into  the  revised  model.  The  two  E-layer  parameters  are  the  critical 
frequency  of  the  auroral  E layer,  ^ . and  the  height  of  the  combined 
solar  and  auroral  E layers,  HE.  The  four  valley  parameters  are  explained 
above.  These  parameters,  with  the  exception  of  two  cases  of  HE,  show  a 
higher  correlation  with  K than  with  K^.  All  have  slightly  lower  a's, 
indicating  that  a specification  with  K will  result  in  better  predictions 
of  the  parameter. 

The  intercepts  and  slopes  of  the  least-squares  linear  fits  for  the 
six  parameters  that  were  used  in  the  model  modifications  are  given  in 
Table  4-3.  In  the  modified  model  these  linear  fits  are  assigned  to  the 
average  x ot  each  solar-zenith  angle  group.  Interpolation  is  used  to 
obtain  the  parameters  for  other  values  of  solar-zenith  angle  (see 
Appendix  C) . 

The  first  parameter  used  in  the  revised  model  is  the  critical 
frequency  of  the  auroral  E-layer.  The  RADC  code  predicts  the  auroral 
E-layer  and  the  solar  E-layer  separately  and  combines  them  afterward 
in  subroutine  ECOMB.  The  critical  frequency  of  the  solar  E predicted 
in  the  RADC  model  was  subtracted  from  the  data,  and  predictions  were 
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CORRELATION  OF  IONOSPHERIC  PARAMETERS  WITH  K FOR  FIVE  RANGES  OF  SOLAR-ZENITH  ANGLE, 


CORRELATION  OF  IONOSPHERIC  PaR-AMETERS  WITH  K FOR  FIVE  RANGES  OF  SOLAR-ZENITH  ANGLE 


SLOPES  AND  INTERCEPTS  FOR  LINEAR  FITS  TO  IONOSPHERIC 
PARAMETERS  FOR  FIVE  RANGES  OF  SOLAR  ZENITH  ANGLE,  X 
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made  on  the  residual,  using  the  method  in  ECOIIB.  Figure  4-3  shows  the 
results  of  linear  fits  of  f E.  versus  K for  various  solar-zenith-angle 
ranges.  The  curves  are  labeled  with  the  average  solar  zenith  angle  for 
each  group.  We  see  that  magnetic  dependence  on  K Increases  with  solar 
zenith  angle.  For  x 60°,  is  nearly  flat;  for  60°  < x 80°,  the 

curve  steepens;  and  for  x ^ 80°  there  is  a very  strong  dependence  on  K. 

The  second  parameter  that  we  predict  is  the  height,  HE,  of  the 
combined  solar  and  auroral-E  layers.  We  first  tried  to  predict  the 
height  of  the  auroral  component  alone  by  subtracting  out  the  predicted 

I solar  layer  height.  This  was  unsatisfactory  because  the  solar-E  height 

prediction  was  unreasonably  high  in  some  cases,  ranging  up  to  160  km. 

t 

j Reasonable  values  were  obtained  for  HE,  but  they  were  not  highly 

correlated  with  K (see  Table  4-2).  As  shown  in  Figure  4-4,  only  a slight 
altitude  variation  was  found  in  the  parametric  analysis.  We  see  that  the 
height  of  the  E layer  generally  decreases  slightly  with  increasing  K and 
increasing  solar  zenith  angle,  except  for  the  solar-zenith-angle  group 
90°  < X ^ 100° , which  shows  the  opposite  trend. 

* i 

t 

All  of  the  valley  parameters  are  predicted  in  the  form  of  ratios  I 

scaled  to  the  f^E^  or  HF^,  which  were  presumed  to  be  well  predicted  \ 

in  the  KADC  code.  (An  exception,  evident  in  Figure  5-5,  is  the  occasional 

occurrence  of  unrealistically  low  altitudes  of  HF^) . We  chose  to  relate 

the  valley  parameters  to  f F„  rather  than  to  f E because  the  correlation 

o 2 o 

coefficient  with  K was  higher.  In  addition,  because  the  Chatanika  data 
were  not  temperature-corrected,  accuracy  is  Improved  if  valley  ratios  are 
scaled  to  the  F region  rather  than  the  E region.  This  effect  is  discussed 
in  Appendix  A.  Figure  4-5  shows  the  scatter  plots  of  ratios  for  the 
critical  frequency  of  valley  one  with  the  best-fit  lines  plotted  on  the 
graph  for  the  two  largest  solar-zenlth-angle  groups.  Figure  4-6  shows 
the  corresponding  scatter  plots  for  valley  two. 

Fig  ire  4-7  shows  depth  of  valley  one  (relative  to  f F„)  as  a 

o 2 

function  of  K for  various  solar  zenith  angles.  Here  there  appears  to  be 
a definite  trend  for  the  deepest  valleys  to  occur  at  low  K values  and 
large  solar  zenith  angles.  Apparently,  for  high  K values  the  F^  region 


58 


M 


FIGURE  4-3  CRITICAL  FREQUENCY  OF  THE  AURORAL  E-LAYER,  foE^,  AS  A FUNCTION 
OF  K FOR  VARIOUS  SOLAR  ZENITH  ANGLES,  x 

is  filled  in  by  particle -produced  topside  E.  The  depth  of  valley  two  as  a 
function  of  K for  various  solar  zenith  angles  is  shown  in  Figure  4-8. 

Here  a similar  trend  occurs,  but  the  x dependence  is  not  as  definite,  and 
for  X 60°  there  is  no  dependence  on  K. 

Figure  4-9  shows  the  height  parameter  of  valley  one  as  a function  of 
K for  various  solar  zenith  angles,  x>  curves  show  a tendency  for  the 

/alley  to  Increase  in  altitude  as  K increases  (x  < 60°),  as  one  might 
expect  from  an  Increase  in  precipitation.  The  only  exception  is  the  line 
for  X 60° , which  is  based  on  a small  data  set  and  was  not  statistically 
slgnif leant. 

Figure  4-10  shows  the  height  ratio  of  valley  two  as  a function  of  K 
for  various  solar  zenith  angles,  x*  These  curves  seem  to  have  a somewhat 
consistent  trend,  with  the  height  of  the  valley  decreasing  as  a function 
of  K for  small  solar  zenith  angles  and  gradually  switching  to  an 
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FIGURE  4-4  HEIGHT  OF  THE  COMBINED  SOLAR  AND  AURORAL  E-LAYERS  AS  A 
FUNCTION  OF  K FOR  VARIOUS  SOLAR  ZENITH  ANGLES,  X 

increasing  height  of  the  valley  as  a function  of  K for  large  solar 
zenith  angles. 

4.3  Profile  Shapes 

The  version  of  tho  RADC-POLAR  model  that  is  the  basis  for  the 
current  modification  provided  four  types  of  profiles.  These  four  types 
were  generated  from  the  eight  ionospheric  parameters  that  were  predicted 
as  functions  of  the  particular  time,  location,  and  magnetic  activity. 

The  elgh*;  parameters  are  height  and  maximvun  electron  density  of  the  E, 
and  layers,  and  semi-thickness  of  the  E and  layers. 

In  order  to  simulate  the  electron  density  profiles  that  are  found 
in  the  auroral  region  we  found  it  necessary  to  add  four  new  profile 
' types  to  the  four  that  existed  originally.  The  main  change  in  profile 
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FIGURE  4-7  RATIO  OF  MINIMUM  PLASMA  FREQUENCY  AT  VALLEY  ONE  TO  CRITICAL 
FREQUENCY  OF  F2  LAYER,  foVi/foF2,  AS  A FUNCTION  OF  K FOR 
VARIOUS  SOLAR  ZENITH  ANGLES,  x 


FIGURE  4-8  RATIO  OF  PLASMA  FREQUENCY  AT  VALLEY-TWO  MINIMUMS  TO  CRITICAL 
FREQUENCY  OF  F2  LAYER,  foV2/foF2-  A FUNCTION  OF  K FOR 
VARIOUS  SOLAR  ZENITH  ANGLES,  x 
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FIGURE  4-9  RATIO  OF  VALLEY-ONE  HEIGHT  ABOVE  THE  E-LAYER  TO  THE  DIFFERENCE 
IN  HEIGHT  BETWEEN  E AND  F2  LAYERS  AS  A FUNCTION  OF  K FOR 
VARIOUS  SOLAR  ZENITH  ANGLES,  X 


K 

FIGURE  4-10  RATIO  OF  VALLEY-TWO  HEIGHT  ABOVE  THE  E-LAYER  TO  THE 
DIFFERENCE  IN  HEIGHT  BETWEEN  E AND  F2  LAYERS  AS  A 
FUNCTION  OF  K FOR  VARIOUS  SOLAR  ZENITH  ANGLES,  x 


shapes  a -e  those  having  a dense  auroral  E-layer.  These  have  been 
modified  to  look  more  like  the  profiles  observed  at  Chatanlka.  Two  of 
the  new  profiles  exhibit  valleys  between  the  E and  layer,  or  the  E 
and  F2  layer.  The  other  two  represent  an  exponential -type  electron 
density  distribution  without  a valley  that  occurs  when  an  F^^  layer  is 
present  or  when  only  an  F^  layer  is  present. 

The  four  profile  types  in  the  original  code  were; 

(1)  A Single  F2-Layer  Profile 

This  type  of  profile  is  a standard  type  that  occurs  at 
night  in  midlatitudes.  A parabolic  layer  is  used  to  represent 
it.  Use  of  this  profile  remains  unchanged. 

(2)  A Two-Layer  (E  and  F2)  Profile 

When  f^E  is  less  than  ^^^2’  profile  is  a fairly 

good  representation  of  a daytime  midlatitude  ionosphere 

without  an  F layer.  However,  when  f E is  greater  than 
i o 

f F this  profile  type  does  not  represent  the  data  well, 
o ^ 

In  this  case,  we  have  substituted  an  auroral  two-layer 

profile  (type  7).  The  RADC  layer-generator  PR0F3  is  still 

used  for  f E < f F 
o o 2 

(3)  A Two-Layer  Profile  With  a Valley 

This  profile  can  be  used  only  when  f^E  is  less  than 
f^F,.  In  many  cases  it  is  a good  representation  of  a profile 
with  a valley.  The  depth  of  the  valley  is  defined  in  the 
RADC  layer-generator  PR0F3  according  to  the  relative  density 
of  the  E and  F^  layers.  To  avoid  conflicting  valley  occurrence 
patterns,  we  omit  this  profile  in  the  region  where  the  depths 
of  the  valleys  are  being  predicted. 

(4)  A Three-Layer  Profile  Without  Valleys 

When  f E < f F < f F„,  this  profile  is  a good  represen- 
o 01  02 

tation  of  a midlatitude  daytime  profile  with  an  F^^  layer.  It 

also  represents  well  the  profile  with  f F,  > f E and  F F,  > f F-, 

olo  olo2 
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a situation  that  occurs  frequently  in  the  auroral  zone.  For 
cases  when  f F,  < f E in  the  auroral  zone  this  profile  is  a 
poor  representation.  We  have  substituted  an  auroral  three- 
layer  profile  (type  8)  for  cases  of  this  type.  For  a type  4 
we  have  added  the  feature  that  if  a valley  occurs  between 
the  E and  F or  between  the  F and  F layers,  it  will  be 

X X w 

represented.  The  RADC  layer-generator  PR0F3  is  then  used  for 
this  profile. 

The  four  auroral  profile  types  that  have  been  added  are; 

(5)  An  Auroral  Two-Layer  Profile  With  a Valley 

This  profile  occurs  when  the  solar  E is  less  dense 
than  the  F^  layer  but  a particle -produced  auroral  E layer 
is  present.  The  shape  of  this  profile  above  the  F maximum 
is  generally  exponential  up  to  the  valley  minimum  in  these 
cases.  The  new  profile  routine  uses  a combination  of 
exponential  for  the  E layer,  and  Chapman  for  the  F^  layer, 
to  generate  the  profiles,  with  the  valley  depth  and  height 
matched  as  closely  as  possible.  Figure  4-ll(a)  shows  a 
profile  of  this  type.  The  plus  signs  are  actual  data 
points;  the  dashed  line  is  a smoothed  version  of  the  data 
points.  Parameters  found  from  the  data  were  supplied  to 
the  RADC  profile  routine  to  generate  the  dotted-line  profile 
and  also  to  the  new  profile  routine  to  generate  the  solid- 
line profile;  thus  the  SRI  and  RADC  profiles  are  forced  to 
agree  at  the  layer  peaks. 

(6)  An  Auroral  Three-Layer  Profile  With  a Valley  Between 

the  E and  Fj  and  Possibly  Between  the  Fj  and  F2 

This  profile  will  occur  when  the  solar  E is  less  dense 
than  the  F^  layer  and  an  auroral  enhancement  of  the  electron 
density  is  present  in  the  E region.  Figure  4-11 (b)  shows  a 
profile  of  this  type.  There  is  no  pronounced  Fj^  layer,  but 
the  layer  is  visible  with  shallow  valleys  below  and  above 
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the  peak.  A local  maxlmiam  between  the  E and  F layers  is 
referred  to  here  as  an  F^^  layer,  whether  or  not  It  is  a 
conventional  F^^  layer.  A combination  of  exponential  for 
the  E layer  and  Chapman  for  the  F^^  layer  is  used  to  generate 
this  profile,  with  the  valley  depth  and  height  matched  as 
closely  as  possible.  Above  the  layer  the  RADC  layer 
generator  PR0F3  is  used.  If  a valley  occurs  it  is  included. 

(7)  An  Auroral  Two-Layer  Profile  (No  Valley) 

This  profile  occurs  when  f E > f F„.  In  this  case, 

o o ^ 

prei'lpltation  has  filled  in  all  signs  of  a valley  between 
the  E and  F^  layers.  An  exponential  shape  above  the  E 
layer  models  the  usual  shapes  very  well.  Very  often  there 
is  a slight  trace  of  the  F2  layer.  Figure  4-11 (c)  shows  a 
profile  of  this  type.  The  Fg  layer  is  visible,  yet  no 
valley  can  be  seen.  To  simulate  this  layer  we  match  the 
electron  densities  at  a point  one-half  the  semi-thickness 
below  the  F^  maximum,  and  use  a combined  Chapman  for  the 
Fg-layer  and  an  exponential  for  the  E layer, 

(8)  An  Auroral  Three-Layer  Profile  (No  Valley  Between  E and  Fj) 

This  profile  occurs  when  f^E  > f^^F^  and  precipitation 
has  filled  in  all  signs  of  a valley  between  the  E and  F^  layers 
Figure  4-11 (d)  shows  a profile  of  this  type.  An  exponential 
shape  above  the  E layer  models  this  profile  type  quite  well. 

A slight  trace  of  the  F^  layer  is  included  by  matching  the 
electron  density  at  a point  one-third  of  the  way  between  the 
height  of  the  F^^  layer  and  the  height  of  the  E layer  with  a 
combined  Chapman  for  the  F^^  layer  and  exponential  for  the 
E layer.  Above  the  F^  layer  the  RADC  layer-generator  PR0F3  is 
used.  If  a valley  occurs  it  is  Included. 

In  designing  the  model,  we  did  not  predict  the  occurrence  of  valleys 
directly.  Instead,  we  assumed  that  ^ a valley  is  present,  its  parameters 
are  defined  by  the  x K dependence  computed  from  the  data  base  of  all 
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profiles  containing  valleys.  A model  profile  is  constructed  from  the 
RADC  F-region  parameters,  the  predicted  E-region  parameters,  and  the 
valley  parameters.  This  profile  is  automatically  subjected  to  a number 
of  tests  to  determine  whether  such  a valley  is  consistent  with  the  E- 
and  F-region  parameters.  Since  the  potential  existence  of  a valley 
depends  strongly  on  the  height  increment  between  the  E and  F layers,  an 
abnormally  low  predicted  height  of  the  F region  was  responsible  for  the 
elimination  of  valleys  in  several  cases.  (For  example,  notice  that  the 
F-reglon  height  is  less  than  200  km  at  the  right  end  of  the  contour  in 
Figure  5-5) . 


4.4  Latitude  and  Local-Time  Variations  of  Auroral-E 


As  indicated  in  Section  4.2,  the  statistical  parameters  extracted 

from  the  Chatanika  data  that  were  used  to  generate  the  ionospheric 

profiles  were  better  correlated  with  K than  they  were  with  K . However, 

P 

for  several  reasons  we  still  wished  to  retain  K as  the  index  of  geo- 

P 

magnetic  activity.  First,  is  a planetary  index  that  is  readily 
available.  Also,  a predictive  model  based  solely  on  K would  require 
measurements  of  the  local  magnetic  variation  over  the  entire  high- 
latitude  auroral  region.  Finally,  the  present  predictive  models  generally 
use  Kp,  and  we  preferred  to  have  our  model  compatible  with  them. 

In  order  to  retain  K , we  needed  to  establish  a correlation  between 

P 

K and  the  local  magnetic  index  K at  hlgh-latltudes . This  was 
P 

accomplished  by  the  analysis  of  one  year  of  magnetic  data  from  16  high- 

latitude  sites.  To  Illustrate  the  results  we  show  in  Figure  4-12  the 

K variation  with  latitude  and  local  time  for  a K value  of  3.  Diurnal 

P 

and  latitudinal  variations,  including  a feature  reminiscent  of  the 

night-side  auroral  oval,  are  readily  apparent  in  the  figure.  The  method 

used  and  the  results  obtained  for  other  values  of  K are  discussed  in 

P 

detail  in  Appendix  B. 


The  establishment  of  a relationship  between  K and  K allowed  us  to 

P 

retain  as  the  index  of  magnetic  activity.  In  addition,  it  resulted 
in  the  Incorporation  into  the  model  an  auroral  oval  with  both  latitudinal 
and  local  time  variations.  A major  shortcoming  of  the  RADC-POLAR  model 
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is  that  for  a specified  value  of  it  predicts  an  auroral  E that  is 
homogeneous  in  both  latitude  and  local  time  throughout  the  oval.  If  we 
assume  that  electron-density  profiles  at  high  latitudes  have  a correlation 
with  K similar  to  that  found  in  the  Chatanlka  data,  we  can  use  the  K- 

versus-K  distributions  to  extend  our  results  to  the  whole  auroral  zone. 

P 

The  results  obtained  do  represent  a realistic  auroral  morphology.  Thus, 
the  modified  auroral  E layer  is  an  Improvement  because  it  varies  with 
magnetic  activity,  corrected  geomagnetic  time,  corrected  geomagnetic 
latitude,  and  solar  zenith  angle. 

The  location  of  the  auroral  oval  has  been  modified  slightly  from 
that  specified  in  the  RADC -POLAR  model.  The  boundaries  of  the  predicted 
nightslde  auroral  oval  were  compared  with  recent  measurements  by  the 
ISIS-2  auroral  scanning  photometer  (Murphree  and  Anger,  1977).  As  shown 
in  Figure  4-13,  the  equatorward  oval  boundaries  were  found  to  be  in 
reasonably  close  agreement.  However,  the  poleward  boundary  of  the 

auroral  oval  in  the  RADC  model  moves  equatorward  with  Increasing  K . 

P 

Although  the  poleward  oval  boundary  is  not  as  well  defined  in 
experimental  data,  the  observations  do  indicate  a poleward  expansion 
with  Increasing  magnetic  activity  (Akasofu,  1977).  The  poleward 
boundary  of  the  photometric  oval  (Murphree  and  Anger,  1977;  Murphree, 
private  communication,  1977)  and  the  40-keV  electron-trapping  boundary 
(Page  and  Shaw,  1972)  are  Indicated  in  Figure  4-13.  We  decided  to 
retain  the  equatorward  boundary  specified  in  the  RADC  model  and  fix  the 
nightslde  poleward  boundary  at  a corrected  geomagnetic  latitude  of  72° . 

The  general  decrease  of  K with  increasing  latitude  provides  a taper  to  the 
I poleward  edge  of  the  auroral  oval  through  the  corresponding  reduction 

f in  electron  density  of  the  auroral  E layer. 

I On  the  dayside  we  kept  the  present  RADC  specification  of  the 

auroral  oval  (dayside  cusp)  because  that  region  generally  is  outside 
the  Chatanlka  field  of  view.  However,  we  chose  to  Include  a dayside 
auroral  E to  the  south  that  is  specified  by  the  K distribution  obtained 
in  our  analysis.  As  discussed  in  Section  3,  precipitation  is  observed 
at  Chatanlka  south  of  the  dayside  oval.  Data  south  of  the  auroral  oval 
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FIGURE  4-13  BOUNDARIES  OF  THE  AURORAL  OVAL  AT  GEOMAGNETIC  MIDNIGHT 

AS  DEFINED  IN  THE  RADC-POLAR  MODEL  AND  AS  DETERMINED  FROM 
ISIS-2  DATA.  The  dashed  lines  indicate  linear  fits  to  ISIS  observations 
(Murphree  and  Anger.  1977;  Murphree,  private  communication,  1977)  and  the 
solid  lines  indicate  the  RADC  boundaries.  Also  indicated  is  the  40-keV 
electron  trapping  boundary  (ETB)  determined  from  ESRO  1 observations 
(Page  and  Shaw.  1972). 
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on  the  daysido  w.M*e  Included  in  the  Chatanika  data  base;  therefore,  we 
feel  that  our  estimates  for  this  region  are  reasonable.  On  the  dayslde 
we  predict  values  of  auroral  E south  to  51°  corrected  geomagnetic 
latitude  using  the  spatial,  magnetic,  and  solar  parameters  described 
above.  South  of  51°  corrected  geomagnetic  latitude  on  the  dayslde 
we  use  1.5  MHz  for  the  auroral  E layer. 

Unrealistically  high  values  of  f^E^  were  predicted  in  the  RADC-POLAR 
model  (Version  dated  2-3-77)  at  all  latitudes  south  of  the  auroral  oval. 

A constant  value  of  3 MHz  was  assigned  to  the  auroral  E-layer  in  this 
region.  We  introduce!  raoi’e  realistic  values  by  assuming  that  south  of 
the  auroral  oval  on  the  nlghtside  there  will  be  no  auroral  activity;  a 
critical  frequency  of  1.5  MHz  has  been  assigned  to  the  auroral  E-layer 
there . 

The  auroral  E layer  in  the  polar-cap  region  remains  unchanged  because 
we  had  no  data  that  could  be  used  to  modify  the  RADC  predictions. 

d . 5 Comparison  of  Chatanika  Data  with  Original  and  Modified  Profile 

Predictions 

In  this  subsection  we  compare  measured  Chatanika  profiles  with 
predictions  from  both  the  RADC-POLAR  ionosphere  and  the  revised  model. 
Because  the  ionospheric  models  are  statistical  mean  models,  one  cannot 
expect  predictions  for  an  arbitrary  specific  day  to  agree  with  actual 
observations.  The  comparison  is  made  simply  to  illustrate  specific 
examples  of  profiles  generated  by  the  models  and  to  show  similarities 
and  difforences  between  them  and  actual  instantaneous  measurements  of 
the  ionospheric  profiles. 

Chatanika  data  from  17  March  1976  are  shown  in  Figure  4-14,  with 
predictions  from  the  RADC-POLAR  model  and  the  SRI  modification  of  it . 

The  dashed  line  indicates  the  measured  profile,  the  dotted  line  indicates 
the  RADC  profile,  and  the  solid  line  indicates  the  SRI  modification.  It 
should  be  recalled  that  both  models  use  the  same  Fg-region  parameters. 

The  measured  profiles  have  not  been  corrected  for  temperature  effects  (see 
Appendix  A).  However,  estimated  corrections  at  three  altitudes  are 
indicated  in  the  figure. 
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• 'I'he  uncorrected  critical  frequency  of  the  layer 
as  mea  <ured  is  generally  less  than  predicted.  A 
temperature  correction  results  in  a better  agreement. 

• The  predicted  layer  critical  frequency  is  much 
too  large  at  0500  UT.  Apparently  Chatanika  was  in 
the  trough  at  that  time. 

• The  midnight-sector  (1100  UT)  predicted  critical 
frequency  of  the  auroral  E-layer  is  smaller  than 
actually  occurred. 

• At  1400  UT  the  southern  edge  of  the  auroral  oval  was 
predicted  to  be  north  of  Chatanika,  whereas  Chatanika 
was  actually  in  the  auroral  oval  at  this  time. 

• Predicted  daytime  densities  for  the  auroral  E layer 
are  greater  than  observed  (e.g.,  0200,  2000  and 
2300  UT) . 

Improvements  in  the  profiles  are  seen  in: 

• The  shape  of  the  profiles  having  an  f F.  > f E is 

O 2 O 

improved.  Valleys  were  not  predicted  by  the  SRI 
model  and  they  did  not  occur. 

• The  shape  of  the  profiles  having  an  f E > f F„  is 

o o 2 

Improved.  At  0800  UT  and  1100  UT  the  shape  directly 
above  the  E layer,  although  not  perfect,  is  more 
like  the  actual  profile. 
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5.  SAMPLE  RAYTRACINGS  THROUGH  ORIGINAL  AND  MODIFIED  IONOSPHERES 


In  order  to  demonstrate  the  usability  of  our  modifications  to  the 
RADC -POLAR  Ionosphere,  we  computed  ray  trajectories  along  long-range 
propagation  paths  in  both  versions,  Ray  trajectories  were  computed  with 
a three-dimensional  raytracing  program  provided  by  RADC  (Toman  and 
Miller,  ]977a,b).  Trajectories  were  computed  for  ordinary-mode 
electromagnetic  waves.  The  particular  propagation  path  chosen  was  one 
that  extended  from  the  western  United  States  toward  Europe,  as  shown 
in  Figure  5-1.  Computations  were  made  for  two  times  on  15  April  for  a 
moderate  level  of  solar  activity.  The  first  time  was  across  the  dayside 
auroral  ionosphere  and  the  second  was  across  the  nlghtside  auroral 
ionosphere . 

5 . 1 Dayside  Propagation  Path 

The  ionospheric  contours  of  local  plasma  frequency  derived  from 
the  RADC-POLAR  model  are  shown  in  Figure  5-2.  The  model  was  sampled 
nt  distance  increments  of  400  km  along  the  surface  of  the  earth  for 
the  time  1800  UT,  a sunspot  number  of  50,  and  a of  2.  A conspicuous 
feature  in  the  ionospheric  contour  is  the  dayside  auroral  E layer 
centered  at  a distance  of  3500  km  with  an  enhanced  F region  above  it. 

Ray  trajectories  wer<'  computed  for  a transmitted  frequency  of  20  MHz 
and  initial  elevation  angles  of  0°  to  12°,  in  2°  increments.  Trajec- 
tories are  shown  only  for  a single  hop  to  the  earth's  surface;  no 
ducting  occurred  on  subsequent  hops.  It  can  be  seen  that  the  wave 
launched  at  a 10°  elevation  angle  is  ducted  to  great  distance. 

The  ionospheric  contours  generated  by  the  modified  ionospheric 
model  for  the  same  time  and  solar/geophysical  conditions  are  shown  in 
Figure  5-3,  The  F-reglon  contours  are  similar  except  for  the  addition 
of  a deep  valley  at  a distance  of  about  6400  km.  The  E-region  critical 
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frequencies  are  slightly  smaller,  and  the  feature  at  3500  km  distance  is 
absent.  Surprisingly,  these  differences  result  in  the  propagation  of 
most  of  the  rays  to  greater  distances. 

5.2  Nlghtside  Propagation  Path 

In  order  to  illustrate  the  transauroral  propagation  in  the  evening 

sector,  rays  were  traced  at  0600  UT  on  the  same  date  for  two  values  of  K . 

P 

The  location  of  the  great-circle  path  and  auroral  oval  are  indicated 
in  the  plot  of  geomagnetic  local  time  vs  latitude  shown  in  Figure  5-4. 

Ionospheric  contours  for  a K of  2 are  shown  in  Figure  5-5  and  5-6. 

In  the  original  ionospheric  model  all  rays  are  reflected.  In  the  modified 
model  the  rays  have  apogees  at  higher  altitudes  because  of  the  reduced 
value  for  the  nighttime  midlatitude  auroral  E.  The  lower-elevation  rays 
then  reflect  off  the  topside  auroral  E layer  before  penetrating  over  the 
polar  cap.  The  ray  that  ends  abruptly  at  ~ 6000  km  was  deviated  outside 
the  azimuthal  sector  modeled  (2°  on  each  side  of  the  great-circle  of  the 
great-circle  plane) . 

The  ray  trajectories  for  the  same  conditions  except  for  a K of  4 

P 

are  shown  in  Figures  5-7  and  5-8.  The  increased  K has  resulted  in  a 

P 

wider  and  more  dense  nighttime  auroral  E layer.  Once  again  the  rays  in 
the  original  ionospheric  model  are  simply  reflected.  In  the  modified 
ionospheric  model  the  rays  launched  at  higher  elevation  angles  are 
ducted  before  penetrating. 
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FIGUHE  5-4  PROPAGATION  PATH  IN  FIGURE  5-1  PLOTTED  IN  COORDINATES  OF 
CORRECTED  GEOMAGNETIC  LATITUDE  AND  LOCAL  TIME 
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FIGURE  5-5  IONOSPHERIC  CONTOURS  OF  LOCAL  PLASMA  FREQUENCY  (MHz) 

GENERATED  BY  THE  RADC  MODEL  FOR  THE  SELECTED  PROPAGATION 
PATH  AT  0600  UT,  A SUNSPOT  NUMBER  OF  50,  AND  A Kp  OF  2.  Ray 
trajectories  are  indicated. 


FIGURE  5-6  IONOSPHERIC  CONTOURS  OF  LOCAL  PLASMA  FREQUENCY  (MHz) 

GENERATED  BY  THE  SRI  MODEL  FOR  THE  SELECTED  PROPAGATION 
PATH  AT  0600  UT,  A SUNSPOT  NUMBER  OF  50,  AND  A Kp  OF  2. 

Ray  trajectories  are  indicated. 
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FIGURE  5-7  IONOSPHERIC  CONTOURS  OF  LOCAL  PLASMA  FREQUENCY  (MHz) 

GENERATED  BY  THE  RADC  MODEL  FOR  THE  SELECTED  PROPAGATION 
PATH  AT  0600  UT,  A SUNSPOT  NUMBER  OF  50,  AND  A Kp  OF  4.  Ray 
trajectories  are  indicated. 
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FIGURE  5-8  IONOSPHERIC  CONTOURS  OF  LOCAL  PLASMA  FREQUENCY  (MHz) 

GENERATED  BY  THE  SRI  MODEL  FOR  THE  SELECTED  PROPAGATION 
PATH  AT  0600  UT,  A SUNSPOT  NUMBER  OF  50,  AND  A Kp  OF  4. 

Ray  trajectories  are  indicated. 
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6,  AURORAL  ABSORPTION  MODEL 


As  a separate  and  quite  independent  task,  we  have  developed  a new 
model  of  auroral  absorption  that  has  an  explicit  dependence  on  magnetic 
activity.  Our  objective  was  to  employ  the  same  frame  of  reference,  and 
a comparable  degree  of  averaging,  as  that  used  in  the  model  of  the  £ 
and  F re;;lons  in  order  to  ensure  consistency  between  the  computed  raypaths 
and  the  corresponding  absorption  estimates.  The  rationale  for  develop- 
ment of  such  a model  is  discussed  in  this  section,  together  with  an 
analytic  definition  of  the  model  itself. 

6.1  Background 

The  model  of  the  auroral  ionosphere  described  in  the  preceding 
sections  does  not  include  the  D region,  where  most  of  the  absorption  of 
HF  signals  occurs.  Neither  the  vertical-incidence-ionosonde  technique 
nor  the  Chatanlka  radar  can  provide  the  data  needed  to  develop  such  a 
model,  because  the  electron  densities  in  the  D region  are  normally  too 
small  to  produce  a measurable  signal.  Only  in  raw  instances  of 
extreme  magnetic  activity  (such  as  the  great  magnetic  storm  of  August 
1972)  can  the  radar  measure  D-region  profiles.  Such  data  are  extremely 
valuable  to  aeronomers  as  a source  of  information  on  the  physics  and 
chemistry  of  the  lower  ionosphere,  but  they  are  clearly  not  sufficient 
for  development  of  a synoptic  model.  Although  other  techniques  have 
been  used  to  measure  D-region  profiles  under  more  normal  conditions, 
these  data  are  also  too  sparse  for  general  model  development. 

On  the  other  hand,  the  lack  of  a D-reglon  model  is  not  a serious  I 

deficiency  in  an  ionospheric  model  intended  for  use  in  HF  raytracing. 

Because  of  the  relatively  high  collision  frequency  between  electrons  and 
neutral  particles  at  these  altitudes,  the  presence  of  free  electrons  in 

the  D region  affects  only  the  strength  of  the  signal,  and  not  its  | 

direction  of  propagation.  For  this  reason,  absorption  that  occurs  in  the  ^ 
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D region  is  called  "nondeviative . " The  customary  procedure  in  HF 
raytracing  calculations  is  first  to  determine  the  raypaths  by  using  a 
model  of  the  K and  K regions,  and  then  to  estimate  the  absorption  at  the 
locations  where  the  computed  rays  pass  through  the  D region.  By  this 
procedure,  a model  of  the  spatial  and  temporal  variations  of  absorption 
can  be  incorporated  directly  into  the  raytracing  code,  thereby  eliminating 
the  need  for  a D-region  model. 

In  contrast  to  the  difficulties  Inherent  in  measurement  of  D-region 
profiles,  absorption  is  measured  routinely  at  numerous  locations  by  the 
relatively  simple  riometer  technique.  A riometer  measures  the  variations 
in  cosmic  radio  noise  that  has  penetrated  the  ionosphere.  Since  the 
cosmic  radio  flux  is  constant  over  long  periods  of  time,  any  changes  in 
intensity  from  one  day  to  another  at  the  same  sidereal  time  can  be 
interpre  ed  as  variations  in  ionospheric  absorption.  After  correc  tion 
for  antenna  beamwidth  and  refractive  effects,  the  result  is  a direct 
measure  of  one-way  vertical  absorption  at  the  observing  frequency.  The 
values  can  be  normalized  to  the  frequency  of  interest  through  the 
inverse-frequency-squared  dependence  of  nondeviative  absorption. 

Riometer  measurements  constitute  a large  data  base  from  which  a 
number  oT  synoptic  models  of  absorption  have  been  developed.  These  models 
have  proven  quite  adequate  for  predictions  at  low  and  midlatitudes, 
where  absorption  is  under  solar  control.  Two  additional  types  of  absorp- 
tion occur  at  high  latitudes:  polar-cap  absorption  and  auroral  absorption. 
We  are  concerned  here  only  with  auroral  absorption,  which  is  caused  by  the 
precipitation  of  energetic  electrons  in  and  near  the  auroral  oval. 

A model  of  auroral  absorption  developed  by  Foppiano  (1975)  has 
recently  been  adopted  (after  some  minor  but  significant  modifications) 
by  CCIR  (1976).  This  is  a model  of  monthly  median  absorption  that 
combines  the  effects  of  substorms  occurring  over  an  entire  month  at  a 
given  hour  of  the  day  and  for  a given  level  of  solar  activity  (specified 
by  the  12-month  running-average  sunspot  number).  Although  the  median 
absorption  is  implicitly  related  to  magnetic  activity  through  its 
dependence  on  sunspot  number,  the  model  does  not  have  an  explicit 
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dependence  on  magnetic  activity.  Consequently,  the  large  Increases  in 
auroral  absorption  associated  with  intense  magnetic  activity  tend  to  be 
obscured  in  the  averaging  process.  Conversely,  a few  large  substorms 
may  contaminate  the  quiet-day  predictions.  This  model  is,  therefore, 
not  appropriate  for  use  with  an  HF  raytracing  code  in  which  the  structure 
of  the  refracting  regions  is  defined  as  a function  of  magnetic  activity 
in  a three-hour  time  frame. 


An  auroral  absorption  model  developed  by  Hargreaves  (1966)  has  an 
explicit  dependence  on  K^.  However,  because  Hargreaves  was  seeking  a 
relationship  between  the  solar  wind  and  auroral  absorption  on  a global 
scale,  his  data  have  been  ordered  according  to  universal  time.  That  is, 
data  from  several  time  sectors  have  been  combined,  thereby  obscuring 
diurnal  effects,  which  are  as  important  as  magnetic  effects  for  our 
purposes. 

The  auroral  absorption  model  developed  by  SRI  follows  Foppiano's 
formulation  closely,  while  incorporating  an  average  magnetic  dependence 
derived  from  Hargreaves’  curves.  The  final  calibration  is  made  by 
matching  the  magnitude  of  the  absorption  to  a third  model  given  by 
Easier  (1963).  Because  these  three  models  are  crucial  to  the  development 
of  the  SRI  model,  their  principal  features  are  summarized  below. 

6 . 2 Review  of  Existing  Models 

The  three  auroral  absorption  models  reviewed  here  are:  the  Foppiano 
model,  the  Hargreaves  model,  and  the  Easier  model.  They  are  all  developed 
in  terms  of  median  absorption  at  approximately  30  MHz  and  are  all 
empirical  models  based  on  riometer  data. 

6,2.1  The  Foppiano  Model 

The  description  given  below  of  the  model  developed  by  Foppiano 
is  essen  lally  that  found  in  CCIR  (1976). 

Foppiano  developed  his  model  by  fitting  equations  to  published 
riometer  data  from  27  stations  in  the  Northern  hemisphere  to  indicate 
the  spatial  and  temporal  variations  of  the  auroral  absorption.  Mapping 
is  in  terms  of  a coordinate  system  of  corrected  geomagnetic  latitude 
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cp,  corrected  geomagnetic  longitude,  0,  and  corrected  geomagnetic  local 

time,  T.  The  percentage  probability,  that  the  absorption  measured 

at  30  MH/.  exceeds  1 dB  for  a riometer  with  an  antenna  pointing  at  the 

zenith,  is  given  as  the  sum  of  two  separate  terms,  Q, ^ and  Q,  . Those 

Id  Is 

terms  correspond,  respectively,  to  absorption  contributions  that  are 
assumed  to  arise  from  two  different  sorts  of  incident  auroral  particles, 
referred  to  as  "drizzle"  and  "splash"  (Hartz  and  Brice,  1967) , Since 
the  sources,  energies,  and  Incident  directions  of  these  particles  differ, 
the  two  terms  have  different  spatial  and  temporal  statistics. 

Empirical  equations  based  in  part  on  the  known  particle 
statistics  and  on  riometer  data  give 


Q = Q 
^1  ^Id 


+ Q 


Is 


(6-1) 


with 


and 


Q,.  = 21  d .d„.d„.d.^.d^ 


Id 


cp  T ''R  “0 


(6-2) 


(6-3) 


The  subscripts  cp,  T,  R,  0 and  M relate  to  terms  which  have  a principal 
dependence  on  latitude,  time  of  day,  solar  activity,  longitude  and 
season,  respectively.  The  terms  are  given  in  Eqs.  (6-4)  through  (6-13): 


d = exp 

9 


- (cp  - cp  ) 
2a  ^ 


with  cp  in  degrees 


(6-4) 


where  cp  = 68  (1 


0.0004  R, „),  in  degrees,  and  a = 3 (1  + 0.004  R,_), 

cp  X2 


in  clegre<;s.  ^2  month  xninnlng  mean  sunspot  number. 


. 1 m 

S = — TsTt 


(6-5) 
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where  T = 10  (1  - 0.002  R,„),  hours 
m 12  ’ 

T is  in  hours  in  the  range  0 to  24  such  that: 

t,  = T for  0 < T < (12  + T ) 

m 

t - T - 24  for  (12  + T ) < T < 24 


d«  = 1 • «•»»  “12 


1 0.58  - 0.42  sin  [0.947  (0  + 85°)]  for  0<  0 < 10 

0*16  for  10  < 0 < 80 

0.58  + 0.42  sin  [l.80  (0  - 130°)]  for  80  < 0 < 180 

0.58  - 0.42  sin  [0.947  (0  - 275°)]  for  180  < 0 < 360 

where  0 is  in  degrees  in  the  range  0-360. 


(6-6) 


(6-7) 


d = 1 - 0.3  sin  3.86  6 
m 


(6-8) 


where  6 is  the  mean  solar  declination  for  the  month 


s = exp 
9 


(cp  - cp^)‘ 


where  = 67  (1  - 0.0006  R ) + 0.3  (1  + 0.012  R ) \t\ , deg. 


t = T - 3 for  0 < T < 15 
t = T - 27  for  15  < T < 24  . 


where  t^  = T 


S.J,  = exp 


for  0 < T < 12 


[-■>:] 


t^  = T - 24  for  12  < T < 24  . 


(6-9) 


(6-10) 


Sr  = 1 ^ 0.009  R^2  . 


(6-11) 


"0  = ‘*0  • 


(6-12) 


s = 1 
m 


(6-13) 


Riometer  data  are  found  to  follow  approximately  a fixed 

relationship  between  Q,  and  A , the  median  absorption.  With  Q as  a 

1 m 1 

percentage  and  A in  decibels 
m 


Qj  + 30 


0,23  Q, 


for  Qj  ^ 1 


for  < 1 . 


(6-14) 


6.2.2  The  Hargreaves  Model 

The  data  for  the  Hargreaves  Model  came  from  six  stations  during 
the  period  April  1963  through  March  1964,  as  shown  in  Table  6-1. 


Table  6-1 

RIOMETER  STATIONS  FOR  HARGREAVES’  MODEL 


Station 

1 

Geographic 

Position 

Geooagnetlc 

Latitude 

Invariant 

Latitude 

South  Pole 

90®00'S.  - 

78-6°S 

Froblaher  Bay 

63°48'N,68°42'W 

75-2°N 

74-5° 

Byrd 

eCf  0l'S,ll9°32'w 

70-6°S 

Great  Whale  River 

55°  20 'N,  77°  50 'W 

66-6°N 

67-8° 

Elghta 

75°15'S,77°10'W 

63-9°S 

Bale  St.  Paul 

47°22'N,7(f  33'W 

58- 8°  N 

60.0° 

Rluaeter 

L-Value  Frequency 
(MHz) 


I ! 
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The  median  absorntion  Is  a.'isumed  to  have  a Gaussian  distribution: 


i 


I 


where  the  magnitude,  A^,  the  latitude  of  maximum  absorption,  X^,  and  the 
width  of  the  absorbing  region,  are  found  from  empirical  curves  that 

relate  absorption  to  magnetic  activity  (K  ) and  universal  time. 

P 

Hargreaves  gives  the  following  formula  for  conversion  between 

A and  Q,  (expressed  as  a percentage): 
m 1 

Q,  = 2200  A (6-16) 

1 m 

6,2  3 The  Easier  Model 

The  model  developed  by  Easier  (1963;  private  communication, 

1977)  Is  based  on  rlometer  data  recorded  at  a chain  of  Alaskan  observing 
stations  over  a five-year  period  (1957-1962)  during  a declining  phase  of 
solar  activity — l.e,,  from  sunspot  maximum  to  sunspot  minimum.  Median 
absorption,  A^,  In  dB  Is  determined  from 

Am  = A(season,  hour) 'B (geomagnetic  latitude). 

The  first  coefficient.  A,  is  found  for  a given  season  (l.e.,  winter, 
summer,  and  equinox)  and  local  time  from  a table  of  values  derived  from 
measurements  at  College,  The  assumption  Is  made  that  the  maximum  always 
occurs  at  65°  geomagnetic  latitude  (l.e.,  near  College),  and  that  auroral 
absorption  occurs  only  In  a region  extending  10°  on  either  side  of  the 
peak.  The  second  coefficient,  B,  is  a triangular  function  that 
approximates  the  geomagnetic  latitude  variation. 

The  median  absorption  values  making  up  the  table  of  coefficients, 
A,  were  derived  from  the  College  data  (Easier,  1963).  These  values  have 
had  corrections  applied  that  make  them  directly  applicable  for  use  In 
estimating  one-way  vertical  absorption.  The  data  have  been  corrected  for 
the  effect  of  the  broad-beam  rlometer  antenna  assuming  a uniform  distri- 
bution of  cosmic-noise  intensity  over  the  antenna  pattern,  a horizontally 
uniform  and  stratified  absorbing  layer,  and  a cosine-squared  polar  diagram 
for  the  three-element  Yagl  antenna.  Corrections  were  also  made  to  remove 
the  small  effect  of  spurious  absorption  produced  when  the  cosmic  noise  Is 
reflected  from  the  F layer  and  sporadlc-E  ionization.  These  data  will  be 

used  for  calibrating  A In  our  final  model. 

m 
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The  composite  EK-  for  1957-1962  was  about  21  (average  K of 

P 

2.6).  Easier  suggests  providing  a K dependence  by  multiplying  A by 

p n 

the  ratio  of  rK^/21.  We  have  estimated  that  the  average  sunspot  number 
for  this  period  was  120. 

6.3  The  SRI  Model 

The  formulation  of  the  model  is  accomplished  in  three  distinct 

steps.  First,  we  modify  Fopplano's  equations  that  define  the  spatial 

and  temporal  locations  of  the  absorbing  regions,  converting  their 

dependence  on  sunspot  number  to  a dependence  on  K^;  the  magnitude  of 

the  maximum  absorption  Is  not  altered  in  this  step.  Next,  we  convert 

Q,  to  A , where  A continues  to  represent  the  monthly  median  value 
in  n 

because  the  magnitude  of  remains  unchanged.  A new  conversion  formula 

is  derived  that  fits  the  data  better  than  does  Fopplano's.  Finally,  the 

magnitude  of  the  absorption  is  expressed  as  a function  of  K by  multiply- 

P 

Ing  the  monthly  median,  A , by  a factor  derived  from  Hargreaves'  and 

Rl 

Easier' s data. 


6.3.1  Definition  of  Absorbing  Regions 

The  expression  for  and  [Eqs.  (6-1)  through  (6-3)] 

from  Fopplano's  model  are  Incorporatee  in  the  SRI  model  without  change. 


However,  the  factors  d 


cp-  cp'  T 

to  reflect  a dependence  on  K rather  than  on  R 


and  S.J,  are  redefined  in  the  SRI  model 


We  first  observe 


p 12- 

that  these  factors  are  Gaussian  functions  that  have  been  normalized  to 

their  respective  maxima  so  that  the  modification  does  not  affect  their 
maximum  contributions  to  Q, . 


In  Eq.  (6-4)  (relating  to  d ) , the  latitude,  cp^,  of  maximum 
absorptl  in  and  the  latitudinal  extent,  a , of  the  absorbing  region  for 

cp 

drizzle  precipitation  are  redefined  as  follows: 

67,864  - 0.6256  K 


cp  = 


75.616  - 2.176 


K < 5 
P 


K , K > 5 
P P 


(6-17) 
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3.06  + .276  K , K 
P P 


-.36  + .96  K , K >5 
P P 


(6-18) 


Although  Ha rtz  (1972)  states  that  there  is  no  latitudinal  shift  with 
Kp  In  the  drizzle  zone,  the  above  equations  correspond  to  those  given 
in  Eq.  (6-4). 


For  the  splash  zone,  the  corresponding  modifications  of  s 


[Eq.  (6-9)]  are: 


66.799  - 0.9246  K + (0.318  + 0.0828  K lltl  ; K < 5 
P P ' P 

78.256  - 3.216  K + ( . 288  K - .708) Itl  ; K > 5 
P P P 


(6-19) 


t = T - 3 for  0 < T < 15 

t = T - 27  for  15  < T < 24 

T = corrected  geomagnetic  local  time  (CXjT). 

Maximum  absorption  in  the  splash  zone  occurs  at  2300  CGT.  The  width  of 
the  splash  zone  is  defined  by  a [Eq.  (6-18)],  Identical  to  that  of  the 
drizzle  zone. 

The  only  other  factors  of  Eq.  (6-1)  that  have  been  altered  are 

d,j,  and  s,^,,  the  times  of  maximum  drizzle  and  splash  precipitation.  For 

drizzle,  the  value  of  T in  Eq.  (6-5)  is  changed  to 

m 


9.9  - 0.46  K for  K < 5 
P P 


15.6-1.6  K for  K >5. 

P P 


(6-20) 


For  a Kp  of  4 this  zone  is  centered  around  0800  as  Hartz  (1972)  suggests. 

For  lower  K values  this  zone  peaks  later  and  for  higher  K it  peaks 
P P 

earlier.  The  time  of  maximum  splash  precipitation  has  been  changed  to 

2300  CGT,  as  recommended  by  Hartz.  Equation  (6-10)  is  modified  as 

follows : 


1 


JT*  1 

|t  - 24  11  < T < 24. 


(6-21) 


It  should  be  noticed  that  the  factors  d and  s , which  express 

R R 

a direct  dependence  on  solar  activity,  have  not  been  altered. 


6.3.2  Conversion  from  Qj  to  A„ 

From  the  above  expressions  the  value  of  can  be  calculated. 

The  conversion  from  Q,  to  the  monthly  median  value  of  A remains  to  be 

1 m 

explained.  The  conversion  factors  given  by  Fopplano  and  Hargreaves 

differ  substantially,  as  shown  In  Figure  6-1.  In  addition  to  the  curves 

of  Q versus  A , the  figure  shows  Hargreave's  data  and  four  additional 
1 m 

points  derived  from  B.TSler's  data.  Fopplano *s  values  of  A are  larger 

in 

than  the  data  seem  to  Imply.  Hargreaves'  curve  matches  the  data  quite 

well  up  to  Q = 10,  above  which  we  had  no  data  for  comparison;  however, 

the  maximum  value  that  A can  attain  Is  0.41  dB  (for  Q,  = 100),  which  Is 

m 1 

clearly  not  reasonable  In  view  of  the  definition  of  Q^.  (Hargreaves 
cautions  against  extrapolating  his  expression  beyond  Its  range  of 
applicability.)  Therefore,  we  have  accepted  Hargreaves'  expression  as  the 
better  1 it  to  the  data  for  S 10.  At  that  point,  we  have  shifted 
Fopplano 's  curve  to  match  Hargreaves',  continuing  thereafter  with 
Fopplano 's  shifted  curve. 

The  resulting  conversion  equations  are: 


y2200  / 

1 

1.6057  \ SO  J 


Q < 10 


Q > 10. 


(6-22) 
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FIGURE 


i-1  COMPARISON  OF  METHODS  OF  CONVERTING 
THE  PROBABILITY,  Q,  THAT  ABSORPTION 
MEASURED  AT  30  MHz  EXCEEDS  1 dB,  TO 
MEDIAN  ABSORPTION,  A„ 


6.3.3  The  Magnitude  of  Absorption  as  a Function  of  Kp 


The  value  of  A obtained  by  the  above  procedure  can  still  be 
m 

Interpreted  as  a monthly  median.  Although  Its  value  at  a given  location 
may  be  altered  by  the  changes  listed  In  Section  6.3.1,  Its  maximum 
value — that  Is,  Its  value  at  the  locations  and  times  of  peak  absorption — 
differs  from  the  Fopplano  value  only  through  the  modified  conversion 
equation  given  In  Section  6.3.2. 

The  final  step  Involves  direct  modification  of  A by  multlpllca- 

m 

tlon  with  a function  of  K . We  define  the  multiplying  factor  as 

P 


Y„  =1.8  (0.17  + 0.14  K ). 
Kp  p 


(6-23) 


The  form  of  this  expression  was  derived  from  curves  given  by  Hargreaves 
(1966).  The  scaling  factor,  1,8,  was  obtained  by  matching  the  computed 
values  to  Basler's  data  from  College  for  the  winter  season,  using  values 
of  R, „ = 120  and  K =2.6  to  correspond  to  the  average  conditions  under 

l ^ p 

which  the  data  were  gathered. 


The  final  conversion  from  the  monthly  median,  A , to  the 

m 

3-hour  median  corresponding  to  a given  value  of  K Is 

P 


A . 
m 


(6-24) 


6.4  Sample  Calculations 

Several  sets  of  calculations  have  been  made  to  compare  the  new  model 
with  other  models.  Figures  6-2  through  6-4  show  contours  of  median 
absorption  at  30  MHz  on  a corrected-geomagnetlc-latltude/corrected- 
geomagnetlc-tlme  grid.  Input  parameters  are  R, o = 120,  K = 2,6,  the 

Iz  p 

month  of  December,  and  the  corrected  geomagnetic  longitude  of  College. 
Figure  6-2  shows  the  new  SRI  abosrptlon  model.  The  drizzle  and  splash 
peaks  are  both  visible,  having  values  of  0.59  and  0 27,  respectively. 

The  Basler  model  Is  shown  In  Figure  6-3.  The  peaks  In  this  case  are 
0.6  and  0.3.  In  the  Fopplano  model  (Figure  6-4),  the  peaks  are  0.95 
and  0.4.  Fopplano's  model  also  shows  higher  absorption  In  the  polar- 
cap  region  In  the  afternoon  sector. 
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GEOMAGNEt>C 


FIGURE  6-2  CONTOURS  OF  MEDIAN  ONE-WAY  VERTICAL  ABSORPTION  IN  dB 
AT  30  MHz  ON  A PLOT  OF  CORRECTED  GEOMAGNETIC  LATITUDE 
vs  CORRECTED  GEOMAGNETIC  TIME,  USING  THE  SRI  ABSORPTION 
MODEL.  Corrected  geomagnetic  longitude,  260°;  sunspot  number,  120; 
month,  12;  Kn  = 2.6. 


FIGURE  6-3  CONTOURS  OF  MEDIAN  ONE-WAY  VERTICAL  ABSORPTION  IN  dB 

AT  30  MHz  ON  A CORRECTED  GEOMAGNETIC  LATITUDE-CORRECTED 
GEOMAGNETIC  TIME  GRID,  USING  THE  BASLER  MODEL.  Corrected 
geomagnetic  longitude,  260°;  sunspot  number,  120;  month,  12;  Kp  = 2.6. 
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Figures  6-5  and  6-6  illustrate  the  effect  of  a change  in  the  K 

P 

index  in  the  SRI  model.  In  Figure  6-5,  K = 1;  in  Figure  6-6,  K =4, 

P p 

Neither  Basler's  nor  Fopplano's  contours  would  be  affected  by  the  change 

in  K . 

P 

The  next  three  figures  (6-7  through  6-9)  Illustrate  the  potential 
use  of  the  model  in  conjunction  with  raytracing.  The  figures  show 
contours  of  median  30-MHz  absorption  on  a corrected-geomagnetlc-latitude/ 
corrected-geomagnetlc-longltude  grid.  The  three  cases  represent 
absorption  maps  corresponding  to  the  ionospheres  used  in  the  raytracing 
examples  in  Section  5.  The  model  is  designed  for  Incorporation  into  the 
RADC  raytracing  code  if  desired.  A listing  of  the  absorption  subroutine 
is  Included  In  the  computer  listing  of  the  Chatanlka  model  in  Appendix  C. 


FIGURE  6-5 


■'^ECTEo  geomagnetic 


CONTOURS  OF  MEDIAN  ONE-WAY  VERTICAL  ABSORPTION  IN  dB 
AT  30  MHz  ON  A CORRECTED  GEOMAGNETIC  LATITUDE-CORRECTED 


GEOMAGNETIC  TIME  GRID,  USING  A LOW  Kn 


(Kp  = 1). 


Corrected 


geomagnetic  longitude,  260  , sunspot  number  120;  month,  12;  SRI  absorption 


geomagnetic  LU" 


FIGURE  6-7  CONTOURS  OF  MEDIAN  ONE-WAY  VERTICAL  ABSORPTION  IN  dB 

AT  30  MHz  ON  A CORRECTED  GEOMAGNETIC  LATITUDE-CORRECTED 
GEOMAGNETIC  LONGITUDE  GRID  fOR  THE  DAYSIDE  PATH,  WITH 
Ko  = 2.  Sunspot  number,  50;  month,  4;  1800  UT. 


7.  CONCLUSIONS 


The  basic  objectives  of  this  study  were  to  examine  the  characteristic 
behavior  of  electron  density  profiles  measured  by  the  Chatanika  radar  and 
to  incorporate  these  observations  into  the  RADC-POLAR  model.  These 
objectives  were  met,  with  the  limitations  Imposed  by  available  time  and 
funding.  The  purpose  of  this  section  is  to  describe  the  major  conclu- 
sions of  our  work,  as  well  as  some  of  its  limitations,  since  they  may 
be  useful  to  other  investigations. 

The  examination  of  the  Chatanika  data  revealed  that  the  major 
features  of  the  high-latitude  ionosphere  could  be  divided  into  two 
components:  a deterministic  component  due  to  solar  photolonlzatlon,  and 
a variable  component,  presumably  due  to  particle  precipitation.  A 
persistent  feature  in  the  auroral  E region  is  a latitudinally  extended 
belt  of  E-layer  ionization  that  is  associated  with  the  diffuse  visual 
aurora.  Another  characteristic  feature  is  the  absence  of  an  F layer  in 
the  dark  ionosphere.  An  enhanced  F layer  is  occasionally  detected, 
although  it  is  transient  in  nature.  On  the  dayslde  the  ionosphere  is 
solar-controlled,  but  enhanced  E layers  are  observed  south  of  the  dayslde 
auroral  oval.  This  is  the  location  of  "drizzle"  precipitation  associated 
with  energetic  particle  precipitation. 

The  presence  of  deep  valleys  in  hlgh-latltude  electron  density 
profiles  is  expected  in  two  circumstances.  The  first  is  in  the  dayslde 
"drizzle"  region  and  the  second  is  the  night  sector  during  the  summer. 

In  both  cases  there  would  generally  be  a solar-illuminated  F layer  and 
a particle-produced  E layer.  Elevation-scan  observations  of  the  summer 
sunlit  diffuse  aurora  would  be  of  interest.  Unfortunately,  such 

i 

measurements  have  not  yet  been  attempted.  j 

In  our  modifications  of  the  RADC-POLAR  model  we  left  many  of  its  I 

major  formulations  unchanged.  In  particular  we  did  not  alter  its  I 

specification  of  the  F layer.  This  was  because  of  the  much  larger  data  | 

104  1 


base  of  F region  measurements  used  In  the  RADC  model  and  the  use  In  this 
study  of  Chatanika  data  that  were  not  temperature-corrected.  The 
primaiy  interest  of  our  work  was  the  valley  and  E-region  where  temperature 
corrections  are  relatively  unimportant.  If  a more  detailed  study  of  the 
auroral  F region  is  made  with  the  Chatanika  radar,  then  it  would  be 
advisable  to  include  temperature  corrections. 

Another  major  fe.'.ture  of  the  RADC  model  that  we  did  not  change  was 

the  use  of  as  a determinant  of  magnetic  activity.  Our  study  found 

that  local  K is  a better  indicator  of  auroral  E.  However,  there  is  a 

preference  for  K^,  because  of  its  widespread  use  and  general  availability. 

Wo  were  able  to  maintain  K as  the  determinant  by  establishing  a 

P 

relationship  between  and  high-latitude  K by  analysis  of  one  year  of 

geomagnetic  data.  Presumably,  a direct  relationship  between  K and 

P 

auroral  E could  be  e.'-^tablished  by  analysis  of  a similarly  large  set  of 
electron  density  data . 

We  also  reduced  the  value  of  the  midlatitude  nightside  "auroral" 

K from  3.0  MHz  to  1.5  MHz.  The  accuracy  of  this  revision  should  be 
verified.  Also,  the  RADC  specification  of  the  polar  cap  and  dayside 
auroral  oval  (cleft)  warrant  further  study  because  additional  measure- 
ments presumably  have  been  made  there  since  the  original  formulation  of 
the  model.  None  of  those  regions  fell  within  the  scope  of  this  study. 

Because  of  its  altitude  resolution  of  6 to  10  km,  the  Chatanika 
radar  is  unable  to  detect  structures  of  smaller  scale  size  that  may  be 
responsible  for  spor.idlc  E.  Although  sporadic  E is  Important  to  HF 
propagation,  we  did  not  attempt  to  include  it  in  the  model. 

Our  analysis  is  similar  to  that  used  in  the  development  of  the 

RADC-F*OLAR  model  in  that  both  produce  statistical  mean  models.  We  did 
* 

not  model  small-scale  (<  100  km)  structures  or  phenomena  transient  on  a 

time-scale  smaller  than  the  3-hour  K index.  As  a resiilt,  the  features 

P 

in  our  m(Klel  are  characteristic  of  persistent,  widespread  phenomena,  such 
as  the  diffuse  aurora  and  the  dayside  "drizzle"  zones.  Features  such  as 
auroral  arcs,  the  diffuse  aurora  "boundary  feature,"  and  auroral  substorms 
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are  not  Included.  However,  the  use  of  a local  K that  varies  with  latitude 
and  local  time  for  a specified  K did  result  In  an  auroral  Ionosphere  with 
more  structure  than  that  In  the  original  model,  and  presumably  a more 
realistic  one. 

Another  similarity  between  our  model  and  the  R/NDC  model  Is  that  they 
are  both  empirical.  A deterministic  formulation  of  the  auroral  E may  be 
possible  and  Is  worth  pursuing.  It  would  be  necessary  to  determine  the 
relevant  ionospheric,  magnetospherlc , or  Interplanetary  parameter  that 
specifies  the  spatial  characteristics  of  the  auroral  ionosphere.  Because 
many  of  these  parameters  (e.g.,  electric  fields,  and  incoming  auroral 
energy  distribution)  can  be  Inferred  from  the  Chatanlka  measurements,  the 
same  data  set  used  In  this  study  could  form  the  basis  for  the  formulation 
of  a deterministic  model. 

Finally  it  seems  appropriate  to  test  the  validity  of  the  raytracing 
methods  by  conducting  a propagation  experiment  in  which  transauroral 
transmission  is  measured  through  an  ionosphere  observed  simultaneously  by 
Chatanikn  elevation-scan  measurements.  Also,  our  model  had  a large 
spatial  scale  size  (~  100  km),  similar  to  that  used  in  conventional 
ra\  tracing  codes.  It  may  be  enlightening  to  examine  simulated  propagation 
thi ough  :in  ionosphere  specified  to  the  10-kra  resolution  of  the  radar 
measurements  to  ascertain  the  importance  of  discrete  structures. 

The  primary  end  product  of  our  work  is  the  modified  ionospheric 
model.  The  major  components  of  it  are  listed  in  Appendix  C,  and  a 
complete  operating  program  is  being  sent  to  RADC  separate  from  this 
report. 
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Appendix  A 


TEMPERATURE  CORRECTIONS  TO  uHATANIKA  MEASUREMENTS 

As  discussed  In  Section  2.1,  the  Chatanlka  measurements  of  electron 
density  are  normally  processed  under  the  assumption  that  the  radar  cross 
section  Is  a constant.  In  fact,  the  radar  cross  section  depends  on  the 
ratios  of  two  Ionospheric  plasma  parameters.  The  first  Is  the  tempera- 
ture ratio,  T^,  the  ratio  of  electron  temperature  T^  to  Ion  temperature 
T^.  The  second  Is  the  Debye  length  parameter,  a,  a number  linearly 

proportional  to  the  ratio  of  the  ionospheric  Debye  length  to  the  trans- 

2 

mltted  radar  wavelength.  For  the  Chatanlka  radar,  a Is  equal  to 

_3 

13.9  T /N,  where  T is  In  K,  and  N is  in  units  of  cm  . The  electron 
e e 

density,  computed  with  a value  of  unity  assumed  for  T and  with  zero 

2 ^ , 
assumed  for  o , is  referred  to  as  the  uncorrected  density,  N . The 

relation  between  this  density  and  the  true  electron  density,  N,  Is 

N = N'(1  + T + o-^Xl  + 
r 

Because  T is  generally  greater  than  1,  the  uncorrected  electron  density 
r 

Is  generally  an  underestimate  of  the  true  electron  density. 

In  order  to  correct  the  measured  electron  density  profiles,  one 

needs  altitude  profiles  of  the  electron  temperatures  and  ion  temperatures. 

These  temperatures  can  be  measured  by  a detailed  analysis  of  the 

frequency  spectrum  of  the  received  signal.  Because  such  an  analysis  is 

time-consuming,  it  is  not  routinely  performed  and  was  not  done  for  the 

data  in  this  study.  However,  temperatures  have  been  determined  for  other 

studies,  and  typical  profiles  are  shown  in  FigUiC-  A-1.  The  correction 

factors  in  the  E and  F regions  for  these  profiles  are  shown  in  Figure 

A-2.  In  the  E region  (below  120  km)  the  electrons  and  ions  are  Isothermal 

and  the  needed  correction  is  less  than  20%  for  densities  greater  than 
4 3 

5 X 10  /cm  . At  higher  altitudes  the  electrons  are  generally  hotter  than 
the  ions  and  the  correction  factor  Is  larger. 
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TEMPERATURE,  K 

FIGURE  A-1  ALTITUDE  PROFILES  OF  ION  AND  ELECTRON  TEMPERATURES  FOR 
VARIOUS  SOLAR  AND  AURORAL  CONDITIONS.  Data  are  from  Kelly 
(1977). 
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FIGURE  A-2  RATIOS  OF  ACTUAL  AND  UNCORRECTED  ELECTRON  DENSITIES  AT  TWO 
ALTITUDES  FOR  VARIOUS  SOLAR  AND  GEOPHYSICAL  CONDITIONS. 
Correction  factors  were  computed  from  temperature  profiles  in  Figure  A-1. 


A specific  objective  of  this  study  is  to  examine  the  depth  of  the 
valley  (defined  as  f V/f  F ) occurring  between  the  i:  and  F region. 

O O a 

Because  the  correction  factor  is  changing  with  altitude,  the  valley  depth 
that  we  infer  differs  from  the  true  valley  depth.  However,  this  difference 
Is  generally  small,  if  the  valley  depth  Is  computed  relative  to  the  F 
region  as  was  done  in  this  study.  An  estimate  based  on  typical  valley 
and  F-region  densities  and  the  temperature  profiles  shown  in  Figure  A-1 
found  differences  ranging  between  an  11%  underestimate  for  the  case  of  an 
active  night,  to  a 4%  overestimate  In  the  case  of  a quiet  night. 
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Appendix  B 


DIURNAL  VARIATION  OF  HIGH-LATITUDE  K 

FOR  VARIOUS  K RANGES 
P 


As  discussed  in  Section  4.2,  the  Chatanlka  electron-density- 

profile  parameters  were  found  to  be  better  correlated  with  the  local 

magnetic  index,  K (College),  than  with  the  planetary  magnetic  index, 

Kp.  In  order  to  retain  the  use  of  as  the  primary  magnetic  index, 

it  became  necessary  to  establish  a relation  between  it  and  K at  high 

latitudes.  With  such  a relationship  we  were  able  to  compute  an  expected 

\'alue  of  K at  any  desired  corrected  geomagnetic  time  and  corrected 

geomagnetic  latitude  for  a given  K . 

P 

Before  describing  the  correlation  that  we  determined,  it  may  be 

useful  to  review  the  meanings  of  the  K and  K Indices.  The  K index  is  an 

P 

Indication  of  the  general  level  of  local  variation  in  the  geomagnetic 
field  during  a 3-hour  period.  The  K-indices  are  scaled  from  normal  and 
storm  magnetograms.  The  three  component  (H,  D,  Z)  records  at  each 
station  are  examined  for  a 3-hour  interval.  For  each  trace  an  estimate 
is  made  of  the  range  during  the  interval.  The  largest  of  the  ranges  in 
the  H or  D traces  is  taken  as  the  basis  for  K.  The  Z component  is  no 
longer  used.  A quantization  scale  from  0 to  9 with  unit  intervals  is 
used.  The  meaning  of  the  scale  depends  on  the  latitude. 


The  Kp  index  is  a measure  Intended  to  express  the  degree  of  magnetic 
activity  over  the  whole  earth.  For  this  purpose,  corrected  or  standardized 
values  of  K (K^)  are  prepared  from  12  observatories,  northern  and  southern, 
in  moderately  high  geomagnetic  latitudes  (see  Table  B-1).  The  mean  of  the 
K values  from  the  12  stations  is  Kp.  It  is  a 3-hour  index  and  is  quantized 
from  0 to  9 in  intervals  of  one-third.  Because  none  of  the  stations  used 


in  compiling  the  Kp  index  are  in  the  auroral  zone,  Kp  is  primarily  a mid- 
latitude  index.  It  is  affected  by  both  the  Intensity  of  the  ring  current 
and  the  intensity  and  location  of  the  electrojet.  For  this  reason  Kp  is 


Table  B-1 


OBSERVATORIES  WHOSE  DATA  ARE  USED  IN  THE  COMPUTATION  OF  K 

P 


Observatory 

Geomagnetic 
North  Latitude 
(deg) 

Geomagnetic 
East  Longitude 
(deg) 

Comments 

Sitka 

60.0 

275.3 

Meanook 

61.8 

301.0 

Agincourt 

55.1 

347.0 

Up  to  1969 

Ottawa 

57.0 

351.5 

Since  1969 

Fredricksburg 

49.6 

349.8 

Hartland 

54.6 

79.0 

Eskdalemuir 

58.5 

82.9 

Lerwick 

62.5 

88.6 

Witteveen 

54.1 

91.2 

Wlngst 

54.6 

94.1 

Rude  Skov 

55.9 

98.5 

Lovo 

58.1 

105.8 

Since  1954 

Amberly 

-47.7 

252.5 

Toolangi 

-46.7 

220.8 

Since  1972 

not  a Kood  quantitative  indicator  of  the  level  of  auroral  activity. 

Although  it  is  a poor  indicator,  it  does  have  some  relationship  to  high- 

latitude  magnetic  activity.  For  example,  Rostoker  (1972)  claims  that  a 

large  value  of  (2  2+)  is  generally  an  Indicator  of  substorm  activity, 

yet  he  also  shows  that  auroral  substorms  can  occur  during  periods  of 

small  K . 

P 

In  order  to  Investigate  the  relationship  between  K and  K , we 

P 

obtained  individual  K indices  for  high-latitude  observatories  from  the 

World  Data  Center  NOAA . We  also  obtained  indices  for  the  same  period 

of  time.  Data  (both  K and  K ) from  1969  to  1974  are  available  on 

P 
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magnetic  tapes,  whereas  more  recent  data  Is  on  microfilm  and  not  as 
readily  available  for  processing.  Data  from  the  year  1969  were  used  in 
this  study. 

In  1969  there  were  16  stations  north  of  53°  corrected  geomagnetic 
latitude  collecting  data  for  K indices.  Table  B-2  lists  the  stations 
used  with  their  code  name,  geographic  coordinates,  and  corrected  geo- 
magnetic coordinates.  The  best  coverage  for  this  year  was  in  the 
European  Sector;  the  coverage  in  North  America  was  rather  sparse.  The 
latitudinal  coverage  is  sparse  above  70° . 

For  each  station  the  diurnal  variation  in  K was  determined  in 

relation  to  specific  ranges  of  K . The  K ranges  used  were: 

P P 

0<K  <2;2<K  <4;4<K  <6;6<K.  If  the  planetary  index  K 
P P P P P 

fell  within  the  range,  the  corresponding  K for  each  station  were 

accumulated  to  get  a yearly  average  K for  that  station  and  for  each  of 

the  3-hour  periods  in  the  day. 

All  of  the  times  were  adjusted  to  corrected  geomagnetic  time.  The 
assumption  was  then  made  that  the  magnetic  activity  measured  by  K is 
dependent  only  on  corrected  geomagnetic  time  and  corrected  geomagnetic 
latitude.  If  the  assumptions  were  true,  stations  at  the  same  latitude 
would  liave  the  same  diurnal  variation  and  level  for  times  shifted  to 
corrected  geomagnetic  time.  We  had  two  sets  of  stations  that  fell  at 
the  same  latitude- -Col lege  and  Kiruna  at  64.8° , and  Tromso  and  Leirvogur 
at  66.8°.  Comparisons  of  the  average  K's  for  both  pairs  showed  similar 
diurnal  variations  that  differed  by  less  than  10%.  The  data  for  each  of 
these  pairs  were  then  averaged  together.  A simple  latitude-weighted 
(1-2-1)  average  was  used  for  latitude  smoothing.  The  final  K time- 
latitude  surface  was  then  generated  by  linear  interpolations. 

Contour  plots  of  these  surfaces  are  shown  in  Figures  B-1  through  B-5 

for  each  range  of  1,  2,  3,  4,  and  5.  Linear  interpolation  is  used  in 

the  model  for  intermediate  K . In  each  case  the  presence  of  the  nlghtslde 

P 

auroral  oval  is  clearly  evident.  The  change  in  size  and  equatorward  shift 

for  Increasing  K is  also  apparent.  For  comparison,  the  latitudinal 
P 

extent  of  the  RAOC  auroral  oval  is  indicated  by  the  vertical  bar  at  mid- 
night in  each  figure. 
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Table  B-2 


1 


STATIONS  USED  IN  DERIVING  THE  TEMPORAL  RELATIONSHIP  OF  K FOR  A GIVEN  K 

P 


Station  Code 


Cape  Cbelyashln  CC 
Dixon  Island  DI 
Tromso  TR 
Godhavn  GO 
Murmansk  MM 
Klruna  KI 
Sodankyla  SO 


Yakulsk  YA 
Domba  s DO 
Nurmljarvl  NU 
Magadan  MG 
Leningrad  LN 
Sitka  SI 
Newport 


Geographic 

North 

Latitude 

(deg) 

Geographic 

East 

Longitude 

(deg) 

77.72 

104.28 

73.55 

80.57 

69.67 

18.95 

69.23 

306.48 

68.25 

33.08 

67.83 

20.42 

67.37 

26.63 

64.87 

212.17 

64.18 

338.30 

63.02 

129.67 

62.07 

9.12 

60.52 

24.65 

60.12 

151.02 

59.95 

30.70 

57.07 

224.67 

48.26 

243.01 

Corrected 

Geomagnetic 

North 

Latitude 

(deg) 

Corrected 

Geomagnetic 

East 

Longitude 

(deg) 

71.58 

174.19 

68.31 

154.69 

66.75 

104.78 

77.62 

41.58 

64.42 

113.87 

64.80 

104.27 

63.90 

108.47 

64.84 

259.65 

66.82 

69.53 

57.58 

198.85 

60.08 

91.89 

57.14 

102.96 

53.91 

217.36 

56.22 

107.44 

59.89 

275.95 

55.38 

299 . 46 

NE 


GEOMAGNEt«C 

FIGURE  B-1  VARIATION  OF  K WITH  CORRECTED  GEOMAGNETIC  LATITUDE  AND 
CORRECTED  GEOMAGNETIC  LOCAL  TIME  FOR  A Kp  OF  1.  The 
boundaries  of  the  RAOC-POLAR  specification  of  the  auroral  oval  are  also 
irKlicated. 


FIGURE  B-3  VARIATION  OF  K WITH  CORRECTED  GEOMAGNETIC  LATITUDE  AND 


^CTeo  GE0MAGN6T»C  ■T' 

FIGURE  B-4  VARIATION  OF  K WITH  CORRECTED  GEOMAGNETIC  LATITUDE  AND 
CORRECTED  GEOMAGNETIC  LOCAL  TIME  FOR  A Kp  OF  4.  The 
boundaries  of  the  RADC-POLAR  specification  of  the  auroral  oval  are  also 
indicated. 
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Appendix  C 


DFi’AILED  SPECIFICATION  OF  COMPUTER  CODE  FOR  THE  CHATANIKA  MODEL 

The  RADC-POLAR  Ionosphere  generator,  subroutine  SUBPOLR,  has  been 
modified  to  Incorporate  the  results  of  our  study  of  electron  density 
profiles  observed  by  the  Chatanlka  radar.  A flowchart  Is  shown  In 
Figure  C-1.  The  modifications  were  designed  to  leave  as  much  of  the 
program  as  possible  unchanged,  since  a great  deal  of  study  and  data 
analysis  has  gone  Into  Its  development.  As  discussed  In  Section  4.1,  we 
have  modified  the  specification  of  the  auroral  E,  and  have  added  four 
new  valley  parameters.  These  parameters  are  generated  In  subroutine 
CHTNIK. 

In  addition  we  have  modified  the  profile-generating  routine 
PROFILE  to  Incorporate  the  profile  shapes  that  were  observed  In  the 
Chatanlkn  data.  The  new  profile  generator  Is  called  SRIPRO. 

Four  options  have  been  provided.  The  first  Is  ISRIMOD.  If  It  Is 
set  to  1,  SRI's  modification  will  be  applied;  If  It  Is  zero,  the  original 
program  will  be  executed.  The  next  three  options,  INOP,  INOG,  INOT, 
allow  the  user  to  omit  the  printing,  graphing,  or  tape  output  If  they 
are  set  to  1;  otherwise  they  are  Ignored. 

A complete  listing  of  the  modified  subroutines  Is  given  In 
Figure  C-2.  Also  Included  is  the  listing  of  the  absorption  subroutine 
(SRISORP)  described  in  Section  6. 
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FIGURE  C-1  FLOWCHART  OF  THE  SRI  MODIFICATIONS  TO  SUBPOLAR, 
THE  IONOSPHERE  GENERATOR  IN  RADC-POLAR 
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FIGURE  C-2  LISTING  OF  SUBROUTINES  USED  IN  SRI  MODIFICATION  OF  THE 
RADC-POLAR  IONOSPHERE  GENERATOR 
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FIGURE  C-2  (Continued) 
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FIGURE  C-2  (Continued) 
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FIGURE  C-2  (Continued) 
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FIGURE  C-2  (Continued) 
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FIGURE  C-2  (Continued) 
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VERTICAL  ABSORPTION  — MEDIAN 
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FIGURE  C-2  (Concluded) 
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